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ABSTRACT 


Numerous  methods  have  been  proposed  for  testing  the  hypothesis 
of  zero  three-factor  interaction  in  contingency  tables.  These  methods 
are  incorporated  into  a  single  computer  program  which  is  designed  to 
calculate  the  test  statistic  by  any  or  all  of  the  proposed  techniques, 
and  which  also  provides  estimators  of  the  interactions  and  correspond¬ 
ing  two-sided  confidence  intervals.  The  program  is  written  in  FORTRAN 
63  for  the  CDC-3600  computer,  and  will  analyze  interactions  in  contingency 
tables  of  dimensions  2  <  r  <  5,  2  <  s  <  5,  2  <  t  <  16,  for  r  <_  s  <  t  . 


COMPUTER  ANALYSIS  OF  THREE-FACTOR  INTERACTIONS 
IN  CONTINGENCY  TABLES 

Marvin  A.  Kastenbaum  and  Dennis  Kuba 


1.  Introduction 

The  analysis  of  three-factor  interactions  in  contingency  tables  has  been  the 
subject  of  numerous  recent  papers  culminating  with  the  admirably  lucid  treatment 
and  summary  given  by  Goodman  [15] .  This  series  of  papers  begins  with  one  by 
Kastenbaum  and  Lamphiear  [8]  in  which  the  authors  present  an  iterative  technique 
for  solving  (r-1)  (  s-1)  (t-1)  simultaneous  fourth-degree  equations.  These  equa¬ 
tions  result  from  an  extension  Bartlett’ s  [  2]  test  of  zero  interaction  in  a  2X2X2 
table  to  the  general  three-way  (rXsXt)  table  as  presented  by  Roy  and  Kastenbaum 
[6].  Subsequent  authors  [10,  11,  12,  13,  14,  15]  have  proposed  alternative  tech¬ 
niques  of  analysis,  each  progressively  simpler  than  the  preceding  one,  but  all 
testing  the  same  hypothesis.  In  every  case;  the  test  statistic  is  distributed 
asymptotically  as  chi-square  with  ( r-1)  (  s-1)  ( t-1)  degrees  of  freedom. 

It  is  the  purpose  of  this  paper  to  describe  a  computer  program  which  has  been 
designed  to  calculate  the  test  statistic  by  any  or  all  of  the  proposed  techniques. 

This  program  also  provides  the  estimates  of  interaction  and  corresponding  simul¬ 
taneous  confidence  intervals  given  by  Goodman  [  15] .  The  program  was  written 
in  FORTRAN  63  for  the  CDC-3600  computer,  and  will  do  analyses  on  contingency 
tables  of  dimensions  2  <  r  <  5,  2  <  s  <  5,  2  <  t  <  16,  for  r  <  s  <  t  .  [Appendix] . 

2.  Tests  of  Hypotheses 

Let  n. ^  be  the  number  of  observations  in  the  ith  row,  jth  column,  and 
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kth  layer  of  a  three-way  contingency  such  that  E  E  E  njjk  =  n>  t3le  tota^ 
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with  similar  relationships  for  the  n  s  .  The  hypothesis  of  zero  three-factor 
interaction  in  an  rXsXt  contingency  table  is  given  by: 


.  nrsknijk  _  nrst.^iJL  ,  2  n 

Ho*  n.  .  n  ..  '  n,  +  n  (2,1) 

isk  rjk  ist  rjt 

for  all  integers  i,  j ,  k  such  that  1  <  i  <  r,  1  <  j  <  s,  1  <  k  <  t  . 

2. 1.  Kastenbaum  and  Lamphiear  f  81 

For  the  rXsXt  table,  Kastenbaum  and  Lamphiear  generalized  an  iterative 
technique  proposed  by  Norton  [  4]  for  handling  the  simultaneous  fourth-degree 
equations  which  arise  in  the  estimation  process.  Under  the  null  hypothesis  (  2.1) , 
estimates  of  the  parameters  may  be  achieved  by  first  solving  for  all  x.^  in  the 
following  systems  of  equations: 


r-l  s-1 

r-l s-1 

(tVsk- .Z  (Z  V"V~V 

(nrSt-Z  Z  V(n«t-V 

s-1  r-1 

r-\  ri 

s-1  r-l  9 

(nisk+jZxi)k)(nrJk+Zxl,k) 

(nist  +  Z1xiJ,'<nrjt  +  ZV 

(  2.1. 1) 


for  all  positive  integers  i  <  r,  j  <  s,  and  k  <t,  and  subject  to  the  constraints 
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be  the 


t-1 

A  x...  =  -x 
k=i  1]k  ' 


for  all  positive  integers  i  <  r  and  j  <  s  .  Let  v... 

ljt  -  —  ijk 


values  of  x  ^  which  satisfy  (  2. 1. 1)  for  all  integers  1  <  i  <  r,  i  <  j  <  s  , 
and  1  <  k  <  t  .  Then 


x2  =  L  £  Z  vfnV(n^  "  vtiJ 

i=l  j=l  k=l 


ijk'  ijk  'ijk' 


(  2.1.  2) 


is  distributed  asymptotically  as  chi-square  with  ( r— 1)  (  s— 1)  ( t— J)  degrees  of 
freedom. 


The  program  solves  equations  (  2. 1. 1)  by  applying  Newton' s  method  of  func¬ 
tional  iteration  to  a  technique  proposed  by  Norton  [  4] ,  and  computes  the  value  of 
the  test  statistic  (2. 1.  2) .  It  then  prints  out  the  number  of  iterations,  the  number 
of  degrees  of  freedom,  Ahe  value  of  the  test  statistic,  and  the  observed  and  ex¬ 
pected  cell  frequencies,  with  identification  for  all  the  cells  in  the  contingency  table. 
A  sample  print-out  is  given  in  Table  1. 


2.  2.  Darroch  [10] 

An  alternative  solution  to  equations  {  2. 1.1)  is  given  by  Darroch.  This 
method  involves  the  iterative  solution  for  6^,  p^.,  and  ijj_  in  the  following 
(rs  +  rt  +  st)  simultaneous  non-linear  equations: 


n 

•Jk  _ 


n 


=  6..  Yj  <p.  ,  i|». , 
jk/^  ki  ij 


(  2.  2.1) 


n 


i.  k 
n 


^ki.^,  ^ij  5jk 


j=l 


(  2.  2.  2) 


n, .  t 

=  ijj. .  ),  6.,  <p,  , 

'  U  k“1  Jk  ki 


n 


(  2.  2.  3) 


U) 

'ki  -  ^ki 


The  iteration  begins  by  setting  <pu.  =  <p\‘,'  =  n,  /n  ,  and  ip..  =  =  /n. 

1  •  K  •  •  K 


i11 

ij  *ij  ij, 


l. 


-(2) 


in  equation  (2.2.1)  and  solving  for  6  ^  =  6^  ,  then  in  equation  (2.  2.  2)  let 
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ib  =  4/  V  and  6Jt  =  and  solve  for  <b.  .  =  <p\  .  In  equation  (  2.  2.3) ,  let 

Tij  Yij  jk  jk  ki  ki  * 

6  ,  =  6^  and  <p,  .  =  <p\  ^  and  solve  for  ik .  =  ib(  ,2^  .  Then  return  to  equation 
jk  jk  ki  ki  ij  ij 

(2.2.1),  letting  ^  ^  and  4^  =  4^  ,  and  solve  for  6jk  =  •  Con- 

,  ,  ,  ,  ,  «  ,  i3,  ,(m)  ,(m-l)  _(m)  Jm-1)  , 

tinue  the  iteration  in  this  way  until  iJj,  .  =  4^j  1  ®jk  =  ^jk  *  anc* 

tpK  -  <pK  to  five  decimal  places  for  all  positive  integers  i  <  r,  j  <  s,  k  <  t 

KI  KI 

At  this  point  calculate 


y2  V  V  V  fn  .(m)  (m)  (m),2i  (m)  ( m)  (m) 

X  =L  L  lj  t  nijk n^jk  *ki  *ij  ]  /n6jk  *ki  *ij  * 


(  2.  2.  4) 


i=l  j=l  k=l 


Darroch  has  shown  that  this  test  statistic  is  distributed  asymptotically  as  chi- 
square  with  ( r-1)  (  s-1)  (t-1)  degrees  of  freedom.  Moreover  (  2.  2.  4)  and  (  2. 1.  2) 
are  algebraically  identical. 


The  program  solves  equations  (  2.  2. 1) ,  (  2.  2.  2) ,  and  (  2.  2.  3) ,  and  computes 
the  value  of  the  test  statistic  (2.  2.  4) .  The  print-out  is  identical  to  that  of  Table  1, 
except  that  the  number  of  iterations  is  generally  smaller.  In  this  sense,  the 
Darroch  procedure  is  superior  to  the  Kastenbaum-Lamphiear  procedure. 

2.  3.  Plackett  f  111 


Let  V =  V/n. .  k  be  the  conditional  probability  that  an  observation  will 
fall  in  the  ith  row  and  jth  column,  given  that  it  is  in  the  kth  layer.  It  follows 


that  the  null  hypothesis, (  2.1^  may  be  rewritten  as 


for  all  integers  1  <  i  <  r,  1  <  k  <  t, 
a  criterion  suggested  by  Woolf  [  5] , 


ijk  ■  aut  1  3 

where  V  “  9rsk9ijk/9isk9rjk  •  Usin9 
Plackett  presents  a  procedure  fortesting  a 


null  hypothesis  equivalent  to  (  2.  3. 1) ,  namely 


(  2.  3.  2) 
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for  all  integers  1  <  i  <  r,  1  <  j  <  s,  1  <_  k  <  t,  where  r  ^  =  log  A.^  .  This 

procedure  is  based  on  the  fact  that  the  maximum  likelihood  estimator  of  log  A_^ 

is  log  d  ,  where  d. =  n. ,,  n  /n.  n  for  all  integers  1  <  i  <  r,  1  <  j  <  s  , 
ijK  ijk  ijk  rsK  is k  rjK 

1  <_  k  <  t  .  Moreover  the  variance  of  log  d  ^  may  be  estimated  consistently  by 

1111 
u , . .  =  +  +  +  • 

1J  nijk  nrsk  nisk  nrjk 


If  R  and  S  are  two  matrices  of  order  (r-1)  Xr  and  (  s-1)  Xs  respectively 
with  rows  orthogonal  to  each  other  and  to  the  unit  vector,  then  the  direct  product, 
[R*S],  is  a  matrix  with  (  r-1)  (  s-1)  rows  and  rs  columns.  The  elements  of 
each  row  of  this  matrix  provide  the  coefficients  of  a  linear  combination  of  the 
logarithm  of  the  frequencies  in  the  kth  layer  of  the  contingency  table.  More 
specifically,  the  matrices  R  and  S  are  formed  as  follows: 


R=  {p  .}  = 

rai 


s=  v  = 


1  for  1<  i_<  , 

-a  for  i  =  a  +  1  , 

0  for  i  >  a  +  1  ,  where  l<i<r,  l<»<r,  and 

1  for  1  <_  j  <  P  , 
t“P  for  j  =  (3  +  1  , 

0  j  >  (3+  1,  where  1  <  j<  s,  1<|3<s  . 


Then  for  each  positive  integer  k  <  t,  a  column  vector  z  ,  is  generated  from  the 

product  of  the  ( r-1)  (  s-1)  X  rs  matrix  [R*S]  with  the  rs  X  1  vector  {log  n  ..  )  . 

ijk 

The  elements  of  z,  are 

k 


r  s 


z.  a  ~  X  X  P  .  <r  log  n. . 

k»3  H,  a  1  61  11 


i=l  j=l 


ijk 


(2.  3.  3.) 


The  asymptotic  distribution  of  z^  =  {z^^}  is  multivariate  normal,  with  dis¬ 


persion  matrix 


Vk  =  [R*S]  D~  ^  [R*S]'  , 
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where  D  *  is  a  square  matrix  of  order  rs  with  elements  { — " —  }  on  the  diagonal 


ijk 

and  zeros  elsewhere.  The  elements  of  V,  are 

k 


ijk 


r  s 

C0V[W  Zka'|3']  Vi'pi  >j/nijk 


(  2.  3.  4) 


for  all  positive  integers  a,  a1  <  r  and  (3,  p1  <  s  .  If  denotes  the  transpose 
of  z  ,  and  V  1  the  inverse  of  V  ,  then  on  the  hypothesis  of  zero  three-factor 
interaction,  (  2.  3.  2) ,  the  statistic 


(  2.  3.  5) 


is  distributed  asymptotically  as  chi-square  with  (  r-1)  (  s-1)  (t-1)  degrees  of 
freedom. 


This  portion  of  the  program  computes  the  elements  of  all  the  vectors  z^  and 

of  their  associated  dispersion  matrices  from  formulas  (  2.  3.  3)  and  (  2.  3.  4) 

for  all  integers  1  <  k  <  t  .  To  evaluate  the  test  statistic,  (  2.  3.  5) ,  (t  +  1)  square 

matrices  of  order  (r-1)  (s-1)  are  inverted.  The  print-out,  Table  2,  displays  all  the 
vectors  z,  and  their  associated  dispersion  matrices  V,  for  all  integers  l<k<t  ; 

t  2K  1  -1  t  -i  _i  t  -1  * 

and  )SlLk  =)&lZk  Vk  Zk*  P  =^k=lVk  ^  9  )?=1  zk  Vk  9  as  wcl1  as  the  value  of  the  test 
statistic,  (  2.  3.  5) ,  and  the  number  of  degrees  of  freedom. 


2.  4.  Goodman  f  14. 151 

The  approach  suggested  by  Goodman  modifies  Plackett' s  method  for  testing 


the  null  hypothesis  (  2.  3.  2)  by  redefining  the  two  matrices  R  and  S  as  follows: 


R  ■  V  * 


for  i  =  a  , 
for  i  =  r  , 

elsewhere,  where  1<  i  <  r,  1  <  a  <  r,  and 
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LfrL.^1 


=  {<Tp) 


1  for  j  =  p 

-1  for  j  =  s 

0  elsewhere,  where  1  <  j  <  s,  1  <  |3  <  s  . 


As  in  Plackett' s  method,  column  vectors  g^  =  {gka^y  are  generated,  where 


r  s 


k£v(3 


=  X  X  p  .  o-  log  n. . 
rai  pj  ljk 


i=l  3=1 


(  2.  4. 1) 


The  asymptotic  distribution  of  each  g^  is  multivariate  normal  with  dispersion 


matrix  U,  ,  whose  elements  are 

k’ 


r  s 


COV[gkap’  gka'  p'^  j^P^iVi  0  Pj  a  p*  j^nijk 


(  2.  4.  2) 


for  all  positive  integers  a,  a'  <  s 
t 

g  =  X  then  the  statistic 


.  If  Uk!  =  Mk,  0  =  [X  Mk]  \  and 

k=l 


-1 


k=l 


Y  =  X  9k  Mk  gk  "  ^  Q? 

k=l 


(  2.  4.  3) 


is  distributed  asymptotically  as  chi-square  with  (  r— 1)  (  s  — 1)  ( t— 1)  degrees  of 
freedom.  Moreover,  expressions  (  2.  3t5)  and  {  2.  4.  3)  are  identical. 

Except  for  the  new  definitions  of  the  matrices  R  and  S,  the  computational 
procedure  using  either  Plackett' s  or  Goodman1  s  methods  appear  to  be  identical. 
In  both  methods  it  is  necessary  to  invert  (t  +  1)  matrices  of  order  (  r-1)  (  s-1)  . 

However,  Goodman  makes  a  significant  contribution,  at  this  point,  by  demon¬ 
strating  that  only  one  (  r-1)  (  s-1)  matrix  and  t  matrices  of  order  (r-1) ,  for 
r  <  s,  need  to  be  inverted  in  calculating  the  test  statistic  (  2.  4.  3)  .  The  re¬ 
duction  in  computing  time  which  results  from  this  modified  procedure  may  be 
appreciable.  The  calculations  are  carried  out  as  follows: 
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(i)  For  every  integer  k,  l<k<t,  define  the  (r-1)  X(r-l)  matrices 

B<k)  -D«w-=t  {n-Hn-}'  ■ 

B*k)=D  fo  Hn  }'  , 

J  nijk  n.)k  lik  1)k 


where  D  and  D  are  ( r-1) X  (r-1)  diagonal  matrices  with 

nuk 


nisk 


elements  n^  and  n.^  respectively,  and  where  {n^}'  and 
{n,.,  }'  are  the  respective  transposes  of  the  (r-1)  X  1  column  vectors 

ijK 

{n^}  and  {n,^}*  for  integers  1<  i  <  r,  and  l<j<s  . 


s — i 

(ii)  Evaluate  ^  +  ^  bJ  ^  ,  an  (r-l)X(r-l)  matrix,  and  its 

( k)  ^ 

inverse  G'  .  1  <  k  <  t  . 


(iii)  Construct  the  submatrices 


M 


(k) 

jj' 


B(k)  (k)  (k)  (k) 

J  J  j 

_(k)  (k)  (k) 

-  Bj  G  B 


for  j  =  j'  =  1,  2, . . . ,  s-1  , 
f  or  j  *  j 1  . 


( iv)  For  every  k,  1  <_  k  <  t,  form  the  square  matrix  of  order  (r— 1)(  s-1) 

2  ( k) 

in  equation  (  2.  4.  3)  from  the  (s-1)  submatrices  Mj^,  each  of  order 
( r-1)  X  ( r-1)  . 

(  v)  Evaluate  the  test  statistic  using  equation(2.  4.  3)  . 

The  print-out  ( Table  3) ,  for  this  section  of  the  program  displays  all  the 

vectors  g^,  their  associated  dispersion  matrices  ,  and 
t  2  t 

S  =  S  gj^  g^,  Q,  ,  as  well  as  the  value  of  the  test  statistic, 
k — 1  k=l 

(  2.  4.  3) ,  and  the  number  of  degrees  of  freedom. 


-8- 


#636 


2. 5.  Goodman  -OH.:  The  2  X  2  X  t  Contingency  Table 


For  the  three-way  contingency  table  with  r  =  2  rows,  s  =  2  columns,  and 
t  >  2  layers,  Goodman  proposes  three  alternative  test  statistics.  Two  of  these 
are  based  on  an  analysis  of  the  cell  frequencies,  and  the  third,  which  is  a  special 
case  of  the  procedure  discussed  in  section  2.4,  is  based  on  an  analysis  of  the 
log-frequencies. 


( i)  Define  d^ 


nllkn22k  u  1  1  1  1 

nl2kn21k  *  k  nllk  n22k  nl2k  n21k 


,2 

Vk  =  dkUk 


and  w^  =  l/v^  ,  1  <  k  <  t 


Then 


x2  =  Yi  dk  wk  "  t  £  dkWk]  2 ^  wk 
k=l  k=l  K  K  k=l 


(2.  5.1) 


is  distributed  asymptotically  as  chi-square  with  ( t  —  1)  degrees  of  freedom. 
Equation  (  2.  5. 1)  may  be  used  to  test  the  hypothesis  of  zero  three-factor  inter¬ 


action  which,  in  the  2  X  2  X  t  table  may 


—  *  - — 


nill  n221  _  nilk  U22k 
°  n!2in211  ni2k  n21k 


2  <  k  <  t  . 


(  2.  5.  2) 


Equation  (  2.  5.  2)  may  be  rewritten  as 


where 


2  <  k  <  t  , 


Ak  = 


nilk  n22k 
ni2k  n21k 


6 Ilk  6  22k 
012k021k 


(  2.5.  3) 


is  a  measure  of  the  two  factor  interaction  in  the  kth  layer  of  the  table.  The 
maximum  likelihood  estimator  of  is  d^,  and  its  variance  can  be  estimated 
consistently  by  v^  . 

Moreover,  Goodman  points  out  that  this  null  hypothesis  may  be  partitioned  into 
the  following  (t-1)  sub-hypotheses: 
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(  given  that  A  ^  =  A  2) 

( given  that  A^  =  A2=A3) 


(  2.5.  4) 


Ht_1:Ai=At  (given  that  A  ^  =  A  ^  =  . . .  =  A^_^)  . 


Each  of  these  sub-hypotheses  ( 1  <  k  <  t)  can  be  tested  using  a  single  de 
gree  of  freedom  as  follows: 


nllvn22v 

Define  d  = 

v  n12vn21y 


1 <  y <k  , 


^  k  k 

and  corresponding  quantities  u^ ,  v  ,  and  w^  .  Let  d^  =  ^  d  w  Jl  w 
*  “  V=1  v  v  V 

and  vk  =  1/2;  w  .  Then  to  test  , 
y=l  Y 


V  *  “i  Y  Y  “ !  V 

1  <  k  <  t ,  in  equation  (  2.  5.  4)  compute 


Xk  ~  tdk+l  "  dk^  ^Vk+1+Vk^  '  (2.5.5) 

t-1  2  2 

The  sum,  2;  X  is  equal  to  X  given  by  (  2.  5. 1) . 
k=l 

When  is  true,  (2.5.5)  is  distributed  asymptotically  as  chi-square  with 

one  degree  of  freedom.  Also  when  is  true,  each  of  the  (k-1)  statistics 
2 

X,  will  have  asymptotically  independent  chi-square  distributions  each  with  one 
degree  of  freedom. 


<U)  Letbk=^-(  ak  =  bj;Uk,  hk=J~,  l<k<t, 

U  i  <2  ,  1  . 

b  =  “  ,  a  =bu,  h  = —  ,  1  <  v  <  k  , 

yd7  Y  YY  y  a  “  * 

Y  Y 

K  =  Tj  bA /E  K >  and  =  K  * 


Y-l  Y  Y  Y=1  V 


Y=1 
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Here  is  the  maximum  likelihood  estimate  of  l/A^  ,  and  its  variance  can  be 


estimated  consistently  by  a^  .  When  the  null  hypothesis  (  2.  5.  3)  is  true 

z2=  i  bkhk-fi  \ m 2/i  v 


k=l 


k=l 


k  kJ 


(  2.  5.  6) 


k=l 


is  distributed  asymptotically  as  chi-square  with  ( t  —  1)  degrees  of  freedom.  Tests 
of  ( 1  <_  k  <  t)  in  equation  (  2.  5.  4)  are  given  by 


Zk  =  tbk+l  "  bk^  ^ak+l  +  ak  ^ 


(  2.5.  7) 


which,  when  is  true,  is  distributed  asymptotically  as  chi-square  with  one 


t-1 


\  2  2 
degree  of  freedom.  Also  ^  Z  =  Z 


k=l 


( iii)  Let 


9k  =  109  dk’  mk  =  ^Uk  >  1  -  k  -  1  » 

g  =  log  d,  m  =l/u  ,  l<v<k, 
Y  Y  Y  Y  —  ~ 

w  lc  if 

* 

g,.  = 


=  £9Ymv/L  mv»  and  \ 

V=1  Y=1 


m 


Here  g^  is  the  maximum  likelihood  estimate  of  =  log  A  ,  and  its  variance 
can  be  estimated  consistently  by  u^  .  The  null  hypothesis,  (2.5.3),  becomes 


Ho:  ri  ■  rk 


for  2  <  k  <  t 


(  2.  5.  8) 


When  (  2.  5.  8)  is  true 


y2  =  m  ]  2/L  mi 

k=l  k=l  k=l 


(  2.5.  9) 


is  distributed  asymptotically  as  chi-square  with  (t-1)  degrees  of  freedom.  This 
statistic,  (  2.  5.  9) ,  is  the  special  form  which  equation  (  2.  4.  3)  takes  when  r  =  2 
and  s  =  2  .  Tests  of  hypotheses  analogous  to  (2.5.  4)  involving  the  log-frequencies 
are  given  by 
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(  2.  5.10) 


Yk  =  [9k+r9kl2/[uk+i  +  uk]  > 


for  every  integer  k,  1  <  k  <  t,  which,  when  is  true,  is  distributed  asymptotically 


V  2  2 

as  chi-square  with  one  degree  of  freedom.  Also,  Y,  =  Y  .  Moreover,  when 

k=l  K 

2  2  K  1  2 

is  true,  and  are  asymptotically  equivalent  to  . 

This  portion  of  the  program  computes  all  the  test  statistics  given  by  equations 

(2.5.1),  (2.5.5),  (2.5.6),  (2.5.7),  (2.5.9)  and  (2.  5. 10).  The  print  *utf 

displayed  in  Table  4,  gives  values  of  d^,  v^,  b^,  a^,  g^,  and  u^  for  l<k<t  , 

2  2  2 

and  X,  ,  Z,  ,  and  Y,  for  1  <  k  <  t  .  In  addition,  the  program  computes  and 

2  2  2  t_1  2  t_1  2  t"1  2 

displays  the  values  X  ,  Z  ,  Y  ,  X  ,  2j  Z ,  ,  }j  Y,  ,  and  the  number  of 

k=l  k  k=l  k=l 


degrees  of  freedom. 

2.  6.  Kullback.  Kupperman.  and  Ku  f  12.131 

The  procedures  outlined  by  Plackett  [11]  and  Goodman  [  14, 15]  for  testing  the 
hypothesis  of  zero  three-factor  interaction,  (  2. 1) ,  are  primarily  generalizations 
of  a  method  proposed  by  Woolf  [5]  in  which  the  test  criterion  is  based  on  a  logit 
transformation  of  the  data.  This  test  criterion  is  one  alternative  to  the  Pearson 
[1]  chi-square  test  of  goodness  of  fit.  Another  alternative  is  the  likelihood- 
ratio  test  criterion  proposed  by  Wilks  [  3] ,  and  investigated  by  Woolf  [  7]  and 
Kullback,  Kupperman,  and  Ku  [12,13]. 

Kullback,  Kupperman,  and  Ku  resort  to  an  information  theory  approach  and 
define  a  minimum  discrimination  information  statistic  [M.D.I.  S.  ].  This  statistic 
is  distributed  asymptotically  as  chi-square  under  the  null  hypothesis,  and  as 
noncentral  chi-square  under  the  alternative  hypothesis,  with  appropriate  degrees 
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of  freedom  and  noncentrality  parameter.  It  has  additive  properties  in  the  sense 
that  it  can  be  analyzed  into  several  additive  components  for  a  hypothesis  which 
is  equivalent  to  the  combination  of  several  hypotheses  of  interest.  Each  such 
component  of  the  M.D.I.S.  is  itself  an  M.D.l.S. ,  and  is  distributed  asymptotic¬ 
ally  as  chi-square  with  appropriate  degrees  of  freedom.  Moreover,  the  M.  D.  I.  S. 
has  the  convexity  property  which  is  useful  in  finding  other  M.  D.  I.  S.  under  certain 
restrictions  and  groupings.  In  particular,  for  nonnegative  real  numbers  a.  and 
b.,  this  property  yields 

n  n  n  n 

X  a  ln(a/b)  >(XaJ  ln<  X  a/X  •  (2.6.1) 

i=l  i=l  i=l  Xi=l  1 

al  a2  a 

Equality  holds  if  and  only  if  “  =  —  =  .  . .  =  t11 

b.  b_  b 

12  n 

In  the  r  X  s  X  t  table,  the  criterion  for  testing  HQ,  (  2. 1) ,  is  given  by 


r  s  t 

=  2  Z  l  l  ni!k  lo9 

i=l  j=l  k=l  J 


n. n.  n  .  n  . 

ilk  l. ,  ,  i.  . .  k 

nn, ,  n.  ,  n  .. 

ij.  l.  k  .jk 


(  2.  6.  2) 


This  statistic  is  distributed  asymptotically  as  chi-square  with  ( r-1)  (  s-1)  ( t-1) 

degrees  of  freedom.  In  any  specific  example  the  convexity  property  may  not  hold. 

This  may  result  in  a  negative  value  for  (  2.  6.  2)  ,  in  which  case  Kullback,  Kupperman 

and  Ku  [12,  pp.  225-226]  recommend  one  of  the  following  alternative  analyses  of 
the  data: 
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row-layer  interaction  with  column 


=  2.Z  E  E  nijkl°9 
1=1  j=lk=l 


ni,kn.i.  yi.k 

nij.n.  jkni.k 


where 


yi.  k  ”  nij.  n.  jk^n.j. 


column-layer  interaction  with  row 

lijk 


2f  = 2  E  E  E  niit  lo9 

i=l  j=l  k=l 


n...  n.  y 
ilk  i. .  .  lk 

n. .  n  ..  n.  , 
i).  .jk  i.  k 


A 

where  y  =  ),  n. .  n.  ,  /n,  ; 

.jk  Ati 


row-column  interaction  with  layer 

2f  =  2  E  E  E  n  log 

1=1  J=1  k=l  1] 


Vn^k^. 

nij.  n.jkni.k 


L 

where  y.,  =  /.  n.  ,  n  /n  ,  . 

ij.  ^=1  i.k  .  jk'  . .  k 


(  2.6.  3) 


(2.6.4) 


(  2.6.5) 


The  program  computes  the  M.  D.  I.  S.  test  statistics  using  equations  (  2.  6.  2) , 

A 

(  2.  6.  3) ,  (  2.  6.  4)  ,  and  (  2.  6.  5) ,  and  prints  out  the  corresponding  values  of  21  as 
well  as  the  number  of  degrees  of  freedom. 

3.  Confidence  Intervals 

Of  all  the  papers  on  the  subject  of  three-factor  interaction  in  contingency 
tables,  only  Goodman1  s  [15]  deals  with  the  problem  of  estimation.  Indeed,  the 
final  section  of  this  paper  is  devoted  exclusively  to  a  discussion  of  methods  for 
estimating  the  magnitude  of  the  three-factor  interaction,  and  for  obtaining  con¬ 
fidence  limits  for  it. 
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Let 


Aibjckd 


(3.1) 


be  the  measure  of  a  particular  three-factor  interaction  for  all  integers  1  <  i  <b  <  r  , 

1  <  j  <  c  <  s,  and  1  <  k  <  d  <  t  .  Depending  on  the  values  of  b,  c,  and  d  , 

r  s  t 

there  can  be  as  many  as  (  2 )  (  2 )  (  2 )  such  three-factor  interactions  in  an 
rXsXt  contingency  table.  In  particular,  for  b  =  r,  c  =  s,  d  =  t,  there  are 
( r— 1)  ( s-1)  ( t-1)  interactions  of  the  form 

Airjskt  =  rijk  “  Fijt  =  1°g<Aijk^AiJt) 


=  log 


(3.2) 


The  maximum  likelihood  estimator  of 


Aibjckd 


is 


and  its  variance  is  estimated  consistently  by 


(  3.  3) 


^ibjckd  Uibjck  +  Uibjcd  * 


1111 

where  u,,  ,  ,  = -  + -  + -  + - 

ibjck  n...  n,  .  n.  .  n,  .. 

ijk  bck  ick  bjk 


(  3.4) 

(3.5) 


Using  these  definitions,  Goodman  proposes  three  alternative  sets  of  approximate 


two-sided  confidence  intervals  for  A.,  .  . 

lbjckd 


aibjckd  *  XT(  P)Sibjc\a 


(3.6) 
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2 

where  x^.(  P)  is  the  [P  X  100]  th  percentile  of  the  chi-square  distribution  with 

T  =  ( r-1)  (  s-1)  ( t-1)  degrees  of  freedom; 

a.u.  .  .  ±  X,(P)S.U.  .  .  ,  (3.' 

lbjckd  1  lbjckd  ’ 

1+P 

where  x,(P)  is  the  [— —  X100]th  percentile  of  the  standardized  normal  dis- 
tribution;  and  for  a  specified  subset  of  interactions,  , 


a 


(w) 


*  V P)  s 


(w) 


(3.8) 


where,  for  any  (  b,  c,  d) ,  w  =  (  b-1)  (  c-1)  ( d-1)  is  the  number  of  elements  in  the 
subset  (w)  =  ( ibjckd)  ;  and 

W  =  W(  B,  C,  D)  =  E  E  E  (b-1)  (o-l)  (d-1)  , 

be  B  cc  Cdt  D 

where  B,  C,  and  D  are  subsets  of  the  sets  of  integers  I  ,  I  ,  I  such  that 

r  s  l 

BCIf  =  {2, 3,. ...r}  , 

O  C  1  =  (2,  3,  . . .  ,  s  }  , 

o 

D  C  It  =  {2,  3,  —  ,t}  ; 

2W-1+P 

and  P)  is  the  ( — 2W~  x  i°°)  percentile  of  the  standardized  normal  dis¬ 
tribution.  It  follows  that  1  <  w  <  W  <  (  r  )  (  a)  (  •  In  particular,  for  the  sub¬ 

set  (  w)  =  (  ir j s kt )  ,  w  =  W  =  T  =  ( r-1)  (  s-1)  ( t-1) ,  and  the  intervals  (  3.  8) 
become 


airjskt  ±  ^irjskt  ’ 


(3.9) 
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2T  —  1  -f-  P 

where  $  (  P)  is  the  ( - - X  100)th  percentile  of  the  standardized  normal 

X  tl 

distribution. 

In  the  final  paragraph  of  his  paper,  Goodman  [15]  discusses  the  relative 
merits  Gf  the  three  alternative  sets  of  confidence  intervals  (  3.  6) ,  (3.  7),  and 
(  3.8).  His  remarks  may  be  summarized  as  follows: 

(  i)  If  any  of  the  intervals  (3.6)  do  not  include  zero,  the  null  hypothesis, 
(2.3.2),  will  be  rejected  at  a  significance  level  ( 1  -  P) ,  when  tested  using 
(2.4.  3). 

(ii)  When  a  specific  set  of  W  three-factor  interactions  is  of  interest, 

W  —  (  2 )  (  2  )  (  2>  ’  ^  anc*  w^en  usual  values  of  P,  [  P  =  .  95  or.  99]  , 

are  used,  (  3.  8)  will  yield  smaller  confidence  intervals  than  (3.6). 

r  s  t 

(iii)  The  probability  is  (  approximately)  at  least  P  that  aU  the  (  2  )  (  2  )  (  2  ) 

intervals  (  3.6)  include  the  corresponding  true  values.  That  is  to  say,  the 

r  s  t 

probability  is  (  approximately)  at  least  P  that  all  the  )  (  2  )  (  ^ )  confidence 
statements  associated  with  (  3.6)  are  correct. 

(  iv)  The  length  of  the  confidence  interval  will  be  reduced  if  (  3.  7)  is  used 
in  place  of  (  3.  6) .  However,  the  probability  that  .all  the  confidence  statements 
are  correct  will  also  be  reduced.  The  expected  proportion  of  correct  confidence 
statements  of  type  (  3.  7)  will  be  approximately  P  . 

(  v)  To  insure  that  aU  confidence  statements  are  correct,  with  probability 
(  approximately)  at  least  P,  use  form  (  3.  6) .  If  it  suffices  to  insure  that  the 
expected  proportion  of  correct  confidence  statements  is  P,  then  (  3.  7)  should 
replace  (3.6). 
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( vi)  The  same  consequences  will  result  when  intervals  (  3.7)  are  used  in 
place  of  intervals  (  3.8). 

For  a  prescribed  set  of  integers  (  b,  c,  d) ,  and  for  all  ( i,  j,  k)  in  the  range 

1  <  i  <  b  <  r,  1  <  j  <  c  <  s,  1  .<_  k  <  d  <  t,  this  portion  of  the  program  calculates 

the  estimates  of  interaction  and  their  estimated  variances  using  equations  (  3.  3) 

and  (  3.  4) .  Moreover,  it  evaluates  the  three  alternative  sets  of  approximate 

two-sided  confidence  intervals  given  by  equations  (  3.  6) ,  (  3.  7) ,  and  (  3.  8) ,  for 

2 

values  of  P  =  0.  95  and  0.  99  .  Values  of  Xip(  P)  in  (  3.  6)  are  calculated  using 
the  Fisher-Cornish  [  9]  approximation. 

Finally,  the  program  is  designed  to  transform  the  estimates  of  interaction  and 
the  corresponding  confidence  intervals  to  their  original  scale  by  taking  the  anti¬ 
logarithm  of  each  of  the  computed  .values.  The  results  of  these  calculations  are 
displayed  in  logarithmic  units  in  Table  5. 
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PROGRAM  h V p o I M T 

CuMi'iCh  /  ALL/u3o(40u)  ♦  I  DEN T  (  10  )  »  IR»J5»<T  »NDEGREES»  INPjT  » OU  T  P 1  T 
CO.V, MON/FLAGS/OPT  I  ON  (  10  )  » LCC-F  L  AG 
TYPE  INTEGER  OPTION 
TYde  INTEGER  OUTPUT 
I  NPUT  =  60 
OUTPUT =61 
10  call  READIN 

IF (OPT  I  ON ( 1 ) .EC.O )  GO  TO  20 
CALL  KASTENBM 

20  IF(OPTION(2)+OPTION(3).EO.O)  GO  TO  30 
CALL  DARROCH 

30  IF (OPTION (A ) • EQ.O )  GO  TO  ^0 
CALL  PLACKETT 

40  IF(0PTI0N(5  ).EC.O)  GO  TO  50 
CALL  goodvani 

50  IF (OPT  ION ( 6 ) .r0.0 )  GO  TO  60 
CALL  GOOOMAN2 

60  IF (OPT  ION ( 7 ) .EQ.O )  GO  TO  70 
C ALu  T aGd  Y  T  aO 

70  I F ( OPT  I  ON ( 8 ) .£0.0)  GO  TO  SO 
CALL  KKK 

80  IF (OPT  ION ( 9  )  .EC.O )  GO  TO  90 

logflag=i 
call  intrvals 

90  IF (OPTION (10). EC.O I  GO  Tn  100 
L  0  G  F  L  A  G  =  0 
CALL  INTRVALb 
TOO  GO  TO  10 
END 


#636 


-21- 


SUBROUTINE  KA5TENBM 

common/ ALL /UBS ( 40 v ) ♦  I  DEN  T ( 1 0  >  »IR»JS»<T ♦  NDEGREES  > INPUT ♦OUTPUT 
DIMENSION  EXPI400) »BIJ( 16) »CIJ( 16) 

TYPE  INTEGER  OUTPUT 
2000  FORMAT ( ]H1 » 10A8  ) 

2 Ou  1  FORMA  M  32HOKAS TENSAUi*i~LAMPri I  E AK  PROCEDURE  ) 

2002  FORMAT (//// ,12h  ITERATION  . 6X *  lShDEGREES  OF  FREED OM*SX* 

*  1 4H  CHI  -UUU  ARE»i5X*2h  R  *  8  X  ♦  2  H  o  »  6  X  ♦  2  ri  I) 

2  003  FORMA  T  (//////*6X»l3.H  CELL  *8X  *  18HU8SERVE  D  FREQUENCY  *  8 X ♦ 

*  1 8HEXPEC  T  ED  FREQUENCY  /»1H  ) 

2  004  FORMA!  (  /  3X  »  I  4  » 18X  *  I  5  *  i  3X  ,  F  1 5  •  6  *  7X  *  3  (  8  X  *  I  2  )  ) 

2005  FORMAT (6X*3(2X*I3)»8X.F15.6»11X,F15.S) 

2006  FORMAT ( 25H0I TERAT ION  COUNT  EXCEEDS  13) 

WRITE  OUTPUT  TADE  OUTPUT *2000* IDENT 
WRITE  OUTDJT  TAPE  OUTPUT *2001 

D£LX=1 .0E-5 
I TERSTOP= 1 00 

nr$t=ir*js*<t 

DO  10  N=1*NR5T 
EXP ( N ) =00  S ( N ) 

10  continue 

I R 1 = I R- 1 
JS1= JS-1 
NX=IR1*JS1*KT 
I TERS=0 
20  NXZER0=0 

ITERS=ITERS+] 

DO  60  I =1 > I Rl 

DO  50  J  =  1  *  JS 1 

SUMC=0.0 

SUmBC  =  0  •  0 

DO  30  <  =  1 *  K.  T 

<PART =  IR* ( Jo*K-JS-l ) 

I  JS*  =  <PART+IR*JS+I 
EXPISK  =  EXP (  I JS<  ) 

IRJK=<PART+IR*J+I7 
EXPRJK=EXP ( IRJK) 

I  RJS<  =  <PART  +  1 P  * JS+ I p 
EXPRS<  =  EXP  (  I PJS* ) 

I JK=KPAPT+IR*J+I 
EXPI JK=EXP ( I J< ) 

Cl  J6  =  1.0/EXPISIA  +  1,0/EXPRJ<+1  .u/EXPRSK+l.U/EXPl  Jn 
Cl J<=1.0/CI JK 

81 JK=EXPI SK*EXPRJ</EXPRS</EXPI JK 
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SUMC  =  SL/',C  +  C  I  JK 
ouMdC=dUM6C+C I JK*6  I J  < 

C I J ( <  >  =C I JK 
d I J ( K ) =B I JK 
30  CONT I  NOE 

h  I  J  =  SUMC  /  bUiXBC 
DO  40  <  =  1  ♦  K T 
KPART=I»*-(JS*<-JS-1  ) 
a  I  JK=CI  J(K)*(  1.0-HI  J*dI  J(K)  ) 

I  F  (  X  I  JK.Lf.  .DcLX)  NXZERO=NXZERC+  1 
IRJSK  =  XP4=>T  +  IC>*JF+I° 

EX  P (  !RJSK)=EXP(IRJ3K)-XIJK 

I JK=<PART+IR*J+I 

EXP (  I JK  )  =EXP (  I JK ) -X  I J< 

I J^K=KPART+IR*J5+I 
EXP(IJSK)=EXP(IJ6K)+XIJ< 

IRJK=KPART+IR*J+IR 
E X P  (  I  R  JK  )  =  EXP  (  IRJO+XI  J< 

40  CONTINUE 
50  CONTINUE 
60  CONT I  NOE 

IF( ITERS. GE.ITERSTOP)  GO  TO  200 
IF (NXZERO.Lr* NX)  GO  Tu  20 
C  H I SQ  =  0  *  0 
DO  7  0  N  =  1  *NRST 
EXPN=EXP ( N ) 

CHlSO  =  CHlSQ-r(OdS  (  N  )  -  EXPN  )  *  (OBS  (  N  )  -EXPN  )  /EXPM 
70  CONT INUE 

write  output  tape  output. 2003 

I  JK  =  0 

rjo  8  0  K  =  1  *  K  T 
DO  80  J  =  1  . JS 
DO  6  0  I = 1  ♦  I R 
I  JK  =  I  J<  +  i 

wRITc  UUTPoT  TAPE  OU I Pu T  ,  2 005  » I ♦ J . < ♦ 03S (  I JK  )  .  E X P (  I J<  ) 
80  CONT I NUE 


WRITE 

UU  T  PuT 

TAPE 

OUTPuT  .2002 

*v  R  i  T  t 

OU  Tru i 

TaPc 

OuTPuT  .2  004.1  TERS.NDEGREES.CHISU.IR. JS.KT 

oO  TO 

3  00 

2  0  j 

/«  R  I  i  E 

uu f  P^ ( 

TAPc 

Ou i Pu T  .200  6* I TERSTOP 

300 

E  Ts  u 
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SUBROUTINE  A3R0CH 

COMMON /ALL, UBSUOu) *  I  DENT (  10  )  »  I  R  *  *>S  » KT »N DEGREES* INPUT  ♦ OUTPUT 
COMMON  /  SUMS /'OB SR  (  60  )  *08SS  (  80  )  *06ST  (  25  )  *OBSRS  (  16  )  »ORSTR  ( 5 )  ♦ 

*  OBSST (*) »0*SP3T 

COMMON /FLAGS/OPT  I  ON  t 10) »  LOGFLAG 

DIMENSION  EXP ( 400 )» THETA (BO)  *PHI (60) *PSI (25) 

TYPE  INTEGER  OPTION 
TYPE  INTEGER  OJIPUT 
2  000  FORMAT ( lnl »  luAo ) 

2001  FORMAT ( lyHODAKROCh  PROCEDURE  ) 

2002  r  ORi'iA  T  (////*  12h  ITERATION  »bX*18nDLGREES  OF  r  RE  EDOM  *  8  X  » 

*  I4h  En I -sOu ARE » 1 3X » 2 n  R»8X»2H  S»6X*2H  I  ) 

2003  FORMAT i ////// »ttX*l3H  CELL  » 6X » 1 HHUftSERV ED  FREQUENCY  * 8X * 

*  18HEXPFCTED  FREQUFNOY  /  •  1  H  ) 

2  004  FORMAT ( /3X* 14* 18 X* 15  *15X»F15.6»7X  *  3 ( 8  X ♦ 12 )  ) 

2005  FORMAT (6X*3(2X,I3 )»8X*F15.6*11X,F15.6) 

2006  FORMAT (25HOlThRAl ION  COUNT  EXCEEDS  13) 

A R I T F  OUTPUT  TAPE  OU TPUT , 2 000  ♦  I  DENT 
WRITE  OUTPUT  TAP-  CirpuT*20Ci 

I TERSTOP=100 

TOLERsU  =  0.0u0(J:3*G.G0U05 

call  getsums 

c 

DO  140  I  =  1  » I R 
OBSST I =0B  SST ( I  ) 

DO  120  J= 1  * JS 
I J= I + I R* ( J- 1 ) 

PSI ( I J)=OBST( I J ) /OBSST I 
120  CONT INUE 

DO  130  <  =  1 » K T 
<  I  =  K  +  K  T  *  (  I  -1  ) 

PHI ( K I ) =ObSS ( K I ) /OBSRS ( < ) 

130  CONT INUE 
140  CONTINUE 
C 

I TERS=0 
200  I STOP  =  1 

I TERS=ITFRS+1 
C 

DO  230  K=1 »KT 
DO  220  J= 1 » JS 
J<*J+JS*  (  < - 1  ) 

SUM  I =0.0 
DO  210  I  =  1  *  I R 
< I =<+KT* ( 1-1 ) 

I J= I + I R* ( J-l ) 

SUM  I = SuM I +PH I  (XI )*PSI ( I J) 

210  continue 

THETA J<  =  OBSR ( J< ) /OBSRS T/ SUV  I 

I F (  (ThETAJK-THETA(JK)  )**2.CE.TOLERSQ)  I  STOP  =  0 
ThETA(J<)=ThETAJK 
220  CONT I NuE 
230  CONTINUE 
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c 


c 


00  2 60  I  =1  *IR 
DO  250  <  =  1  * K  T 
<.  I  =<+<T*  {  I  -1  ) 
oOV j  =  0 • 0 
DO  240  J=1 » JS 
I J=I+IR*( J-l  ) 
j<  =  j+js* ( <- n 

bJMJ=SOMJfPSl (IJ)*THFTA(J<) 

240  CONT I NUE 

PHIM=ubob(i^I  )  /  0  d  o  RST/bu^J 
I F ( (PHI<I-Pnl (<I ) )**2.GE.T0LERSQ)  1ST OP =0 
PHI  (<I  )  =  PH  I <  I 
240  COMT I  N  J  E 
260  CONT I NUE 


DO  290  J  = 1 ♦ J S 
DO  d  6  0  I  =  i  »  I  R 
I J=I  +  IR*( J-l  ) 

S'JvK  =  0.n 
DO  270  <  =  1 » < T 
J<=J+Jb* ( *-l  ) 

<  J  =K  +  KT* ( T - 1  ) 

o U M <  =  SUM K  +  T  ri E  T  A  (  J <  )  *  P H  I  (<I  ) 

270  CONTINUE 

Poll  J  =  OB3  T  l  I  J)  /  OBSf?ST/SUV< 

IF (  (Pol  I J-Pol ( I J)  ) **2 • GE • T GLERbQ )  I3TCP=0 
P 5  I (  I J ) =PS I  I J 
280  CONTINUE 
290  CONTINUE 

IF ( I TERS.GE. I TEPSTOP )  GO  TO  400 
I  F ( I STOP.EO.O )  GO  TO  200 
IF  (OPT  I  ON  (4  )•£<*•  0  )  GO  TO  295 
a-RITE  OUTPUT  TAPE  OuTPuT»2003 
295  CHIoQ=0.0 

DC  300  <  =  1  * < T 
Du  300  J=l»Jo 
JK=J+JS* ( <-l  ) 


DO  300  I  =  1  ♦  IP 
I J=I  +  IR*( J-l  ) 

<  I  =<.  +  < T*  (  I  - i  ) 

EXPIJA=ThETA(jN)*PHI (  n  I )*P5I ( I  J ) *uBSRST 
lJs=I  +  IK*(J-l  +  Jo*(\-l)  ) 

0?,o  I  J<  -=OAS  (  I  JX  ) 

Cril ou  -CHl ou+l ucol J\-EAP I JA ) * (G3SI J<-EXPl JA ) /EXP  I J< 

IF (OPT  I  ON (4 ) .EQ.C )  CO  TO  300 

wkITe  uulPul  TAPE  uuTr'uT»2uJ5»l»J*N*GbblJ\»EXPlJl< 

^no  cont j nuf 

a  R I T  E  OUTpuT  T  APc  OJTPuT»2002 

K  M  L  uuTpuT  TAPl  OuTPuT  *2uU4»lTERu*NDEGREE3»CHlS0»IR»Jb»<T 
GO  TO  500 


400  akite  output  tape  uuTpu t *2U06 * i tfrstop 

rOO  END 
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SUBROUTINE  PLACKETT 

CO  Mil  ON  /AlL/uds(40U)»IDlNT  (1u)»IR»JS»<.T»NDEGR::ES»INPUT  *  OU  TPUT 
DIMENSION  V  (  2  5  6  )  » »/ 1  N  V  (  2  5  6  )  *4(16)  » 2  V  I  N  V  (  1 6  ) 
type  INTEGER  output 

2000  FORMAT ( 1H1 » 10A8 ) 

2001  FORMAT  (20H0PLACK.ETT  PROCEDURE  ) 

2002  FORMAT ( /// /  /  1  y  rl  DEGREls  OF  FREEDOM  *6X,14H  CH I -SQUARE  *  1 7X , 
*  10H  SUM(L*L)  *  1 5 X ♦ 2 H  R»8X»2H  S»RX,2H  T) 

2003  FORMAT (/6X» 15* 15X,F15.6*11X,F15.6»8X,3(8X,I2)) 

2004  FORMAT (3H0<=  12) 

2005  FORMAT I/12H  CONTRASTS  ) 

2006  FORMAT (/20ri  DISPERSION  MATRIX  ) 

2007  FORMAT ( /////lln  Z-VtCTOR  ) 

2008  format ( /  hh  p-matrix  ) 

ftRlTE  OUTPUT  TAPE  OUTPUT ,2000 » IDENT 

WRITE  OUTPUT  TAPE  OUTPUT, 2001 

I R 1 = I R- 1 

JS1=JS-1 

IR1JS1=IR1*JS1 

ZVINVZ=0.0 

DO  30  NA6=1*IR1JS1 

ZVINV(NAB) =0.0 

DO  30  NCD=l»IRiJSl 

NABNCD=NA5+IR1JS1*  (  NCD-1  ) 

V INV (NABNCD)=0.0 
30  CONTINUE 

DO  200  <= 1 , < T 
DO  40  NAR=1 , IP1JS] 

Z ( NAB) =0.0 

DO  40  NCD=1  »  IP  1 JM 

NABNCD=NAb+IRl JSl* ( NCD-1 ) 

V ( NABNCD ) =0.0 
40  CONTINUE 

DO  170  J= 1 , US 
DO  170  I = 1 » I R 
IJK=I  +  iR*(J-l  +  JS-M6-l  )  ) 

OBSI J<=03S(  I  UK ) 

05SL0G  =  L0GF ( CBS  I J< ) 

DO  160  NB  =  1 » JS 1 
S I GMA= 1 • 0 

I F ( J+NB-JS-1 )  60,50,170 
50  S I GMA=-FLOATF ( JS-No ) 

60  DO  180  N A  =  ]  ,  IR1 
PHO=1.0 

IF ( I+NA-IR-1 )  80,70,160 
70  RHO=-FLOATF( IR-NA ) 

80  NAB=IR-NA+IR1*( JS1-NB) 

Z ( NAB ) =Z ( NAo  )+RhC*SI  GMA*05C'|_'*'G 
DO  140  ND=i »JS1 
OMEGA= 1 • 0 
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IF ( J+ND-JS-1 )  100»P0»150 

90  CVEGA  =  -F|_OATF  ( 

100  DO  130  NC  =  l*m 
T  A J= 1  .0 

IF(l-n\C-IR  —  i)  j.20*ilU*l4U 
110  TAU  =  -FLCATF (  I R-NC  ) 

120  NCD=ia~'JC  +  IRl*(  JS1-MD) 

NAoNC L'*NA3+  I  R  1  JS  1*  (  NCD-1  ) 

V  (  NAoi^CD  )  =  9  I  NmdNCD  )  +  R.-.o*o  I  G'*  A*T  Au*Qi-iEG  A /DBS  I  J< 

130  CONTI NU~ 

140  CONTINUE 
150  COMINU1* 

1  5  0  C^mT I  MU- 
170  CONTI  No? 

/.RUE  OuTpn  TAPE  OU 1 PuT , 2 004  ♦  < 
a  RITE  CuTPuF  TAPE  OuiPJT»2005 
CALL  NMPrt  1 1\  T  (  Z  ♦  I  R  1 J  o  1  »  1  ) 

/.  R  I  T  E  OuTPUT  l  A  P  r  OUTPUT *2 005 
CALL  N  i'i  P  R  I  i\  1  (  v  »  I  K  i  J  i  i  »  I  R  1 J  o  1  ) 

C  A  L  u  A  T  I  hV  l  v  *  I  «  Uoi  »v  *0*0*1 rclJSl  ) 

DO  190  NCD=1 >  I  R1JS1 
IM\  NVCD=Z9 I  N 9  <  NCD  ) 

ZCD=Z(NCD) 

DO  180  N A 6  = 1  * IR1JS1 
ZA 5  =  2 ( NAB  ) 

NABNCD  =  \'Ad  +  I  P  1  JS1*  (  NCD-1  ) 

9AriC D=V ( NAoNCD ) 

Z9IN9Z=Z9  INVZ+ZA6*VA5CD*ZCD 
VIN9(i\AdNCDi=\/INV(NAbNCD)+VADCD 
ZVlN9Cu=ZV  INVCu  +  ZAo*'9A5CD 
180  CONTINUE 

c\J  I  -W  (  NCD  )  =4  V  I  NVCu 
190  CONTINUE 
700  CONTINUE 

CALL  ,viATIN9(9IY9*IRlJwl»V*0»D*IRlJS]  ) 

CHI S  P  =  Z  V  I N  V  Z 
DO  210  N C D  =  1  *  I R1JS1 
Z  V  I  N  V  C  D  =  Z  V  I  N  V  (  N  C  D  ) 

DO  210  N A 6  =  1  * m J£1 
NABNCD=NAd+I R 1 JS1*( ncd-i I 

ChIoO  =  uHI  bvj-z  v/  I  NV  (  NA o  )  *91  NV  (  NA5NC0  )  * Z  V  I  NVCD 
710  CONT inue 

/« r i  r e  jo i Pu r  Tape  output  »2UU7 
call  NNiPPINT  (  Z V  I  N V  »  I  R  1 J S  1  *1  ) 
a  R I T  E  OuTPjI  TAPE  0UTPuT*2008 
call  NVPRINT ( V  I NV* IP1 JSl ♦ I  PI Jni  5 
/.rite  output  Tapp  our put, 200? 

AKI  Te  uuTPuT  I  APc  UjIpuI  *200_}*NDEGREEo*CHlbu»ZVliWZ»lR»JF*<.T 
7Q0  END 
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SUBROUTINE  GOODMAN  1 

COMMON / AL L /08S (  40u ) ♦ IUENT (10)*IR*JS*<T  * NDEGREES » I NPUT  *OUTPUT 
DIMENSION  V ( 256 ) »VINV(256)  *Z(16) *  Z V  I NV ( 1 6 ) 

TYPE  INTEGER  OUTPUT 

2000  FORMAT ( lril * 10A8  ) 

2001  FORMAT ( 19HOGOODMAN  PROCEDURE  ) 

2002  FORMAT ( /////19H  DEGREES  OF  FREEDOM  .8X.14H  CHI -SGUARE  *  1  7X  » 
*  10H  SUM ( H*H )  » 1 5  X  »  2  H  R.8X.2H  S.8X.2H  T) 

2003  FORMAT (/6X»I5*15X*F15.6»11X»F15.6»8X»3(8X*I2)) 

2004  FORMAT (3H0<=  12) 

2005  FORMAT ( / 1 2 H  CONTRASTS  ) 

2006  FORMA T ( / 2 UH  DIoPERSluN  MATRIX  ) 

2007  FORMAT (/////11H  G-VECTOR  ) 

2008  FORMAT (/  1 1  n  O-MATRIX  * 

WRITE  OUTPUT  TAPE  OUTPUT ,2 000 ♦ I  DENT 

wRITE  OUTPUT  TAPE  OUrPuT.2001 

I R 1 = I R- 1 

JS1= JS-1 

I  R 1  Js 1  =  I R 1  * JS 1 

ZVINVZ=0.0 

DO  30  NA8  =  1  *  I R 1 JS 1 

Z V  I NV ( N A8 ) =0 • 0 

DO  30  NCD=1*IR1JS1 

NABNCD=NAB+IR1JS1 * ( NCD-1 ) 

VINV(NABNCD)=0.0 
30  CONTINUE 

DO  200  <  =  1 ♦ < T 

DO  40  N A8= 1 » I R 1 J S 1 

Z(NA8)=0.0 

DO  40  NCD  =  1  *  I R 1 JS 1 

NA8NCD=NAB+IR1J31*(NCU-1 ) 

V ( NABNCD ) =0 • 0 
40  CONTINUE 

DO  170  J= 1 ♦  JS 
DO  170  1=1. IR 
I  JK=I  +  IR*(J-1  +  JS*  U-l ) ) 

OBS  I  JK=03S  (  I  Jkr  ) 

OBSLOG=LOGF (03SI J< ) 

DO  160  NB=1 » JSl 
S I GMA=- 1 • 0 

IF(J.EQ.JS)  GO  TO  60 
IF(NB.NE.J)  GO  TO  160 
S I GMA= 1 • 0 

60  DO  150  N A  = 1 » I R 1 
RH0=-1 .0 

IF (  I  . EQ . I R )  GO  TO  80 
IF (NA.NE. I )  GO  TO  150 
RHO= 1 • 0 

80  NAB=NA+IR1*(N5-1 ) 

Z(NAB)=Z(NAb) +RhG*8 I GM A *08 SLOG 
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DO  140  ND=  1  » JS1 
0 VFG4  =  - 1  «  0 

IF(J.EO.JS)  GO  TO  100 

IF1ND.NE.JJ  GO  To  140 
ONEGAs i.o 

100  DO  130  NC  =  1 »  I P  1 
T  Au  =  - 1 .0 

IF ( I ,EQ. IR )  GO  TO  120 
IF  (NC.NE. I  >  GO  TO  130 
T  Au  =  1  •  0 

120  NCD  =  NC  +  IR1*<ND-1  ) 

NA6NCD  =  NAd-*-I  R  1  J S  ]  *( NCD-1  ) 

V  (  iNiAdNC  D  )  =  V  (  NAoi'JC  D  )  -t-RnU*G  1  G  ."T  A  *  T  AO  *  CM  EGA  /  CBS  I  J< 

130  CONTINUE 

140  CONT INUE 
1 50  CONT INUF 
1 f 0  CONT T  NUP 
i 70  CONT I NUC 

n  K  I  I  t  OulFvjT  I  APl  UoTr'oT»2uU4»flv 
aRITE  OUTPUT  TAPE  QjrpoT*2005 
CALL  NMPRINT(Z  ,  IP1JS1  »1  ) 
aRITE  OUTPUT  i APE  OUTpuT*20C6 
CALL  NMPRINT (V* I K 1 JS1 ♦ IRIJdl  ) 

CALL  RiATINVlV*IRlJol*V*(J*D*IRlJi)l  ) 

DO  190  NC D  = 1 ♦ IR1JS1 
ZVINVCD=ZVINV(NCD  ) 

ZCD  =  Z ( NCD  ) 

DO  ] an  M A R  =  ]  *  T  D 1 J P 1 
Z A B  =  Z ( NAB  ) 

NAoNCD=NAb+I R 1 JS] *( NCD-1 ) 

VABCD=V( NAgNCD > 

L\l  I  NVZ=ZV I NVZ+ZAd*VAuCD*ZCD 

V  I NV ( NABNCD ) =  V  I N V ( NAdNCD ) +VAdCU 
Z\/INVCD  =  ZV  INVCD+ZAB*VABCU 

180  CONT INUE 

ZVINV( NCD ) =Z V  I NVC  D 
190  CONT INUE 
200  continue 

CALL  A  T  I  N  V  (vIN\/*IRiJol*V*U*D*  I  R  1 J  S  1  ) 

CHI 50=ZV I MVZ 
DO  210  NCD=1*IP1 JC1 
Z\/INVCD  =  ZVIfW(NCD  ) 

DO  210  NAB  =  !  ♦  1  R i  J F.  1 
NAdNCD=NAd+IK1Uo1*(NCu-1  ) 

Crtiou  =  CnIi>U“Z\/INV(^AD)*\/I^V(NAcNCD)*Z\/INVCD 
210  CONTINUE 

aRITE  OUTPUT  TAPt-  0JTPuT»2007 
CALL  NMPR I  NT (2VINV*IRiJGl*l  ) 
aRITE  OuTPuI  TAPc  u^IhoT»2uU6 
CALL  NRjPk  I  NT  (  \/ 1  imV  *  I  R  1  jG  1  *  I  K  1  JS  1  ) 
aRITE  OUTPUT  TAPE  OUTPUT  *2002 

aRJTe  OUTPUT  TAPE  UuFPUTt200  3»\DFGRE:rG*CHIi:G*ZVI\VZ*TR»Jc»KT 
300  END 
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SUBROUTINE  GOODMAN2 

COMMON /All/UBS(4ow ) * I  DEN  T ( 10) »IR»JS»KT  »NDEGREES* INPUT  *  OU  TPUT 
DIMENSION  P  <  2  56 ) *0(256) *6(80) »C( 16) *G( 16) »PGC 16) 

TYPE  INTEGER  OUTPUT 

2000  FORMAT { 1H1 »10A8 ) 

2001  FORMAT  (  28HOMOD  I  F  I  ED  GOODMAN  PROCEDURE  *////lh  ) 

2002  F0RMAT(/12H  CONTRASTS  ) 

2003  FORMAT  ( /////19H  DtGREnS  OF  FREEDOM  »6X*14H  CHI-SQUARE  *15X  *2H  F 

*  8X»2H  S.8X.2H  T  ) 

2004  FORMAT  (/6X»I5»15X»F15.6*7X»3(8X*  12  )  ) 

WRITE  OUTPUT  TAPE  OU TPUT » 2  GOO  *  I  DENT 
WRITE  OUTPUT  TAPE  OUTPUT»2001 

I R1 = I R-l 

JS1=JS-1 

IR1  IR1  =  IR1*IR1 

IR1 JS1=IR1*JS1 

CHISQ=0,0 

DO  15  NAB  = 1 » I R 1 JS 1 

PG(NAB)=0.0 

DO  10  NCD=1 » I R 1 JS 1 

NA6NCD=NA6+IR1JS1*(NCD-1 ) 

0 ( NA3NCD ) =0  *0 

10  continue 

15  CONTINUE 

DO  200  K  =  1 » K T 

DO  25  N  A6=  1 » IR1JS1 

G(NAB)=0.0 

DO  20  NCD= 1 » I R 1 JS 1 

NABNCD  =  NAB+IR1 JS1*(NCD-1  ) 

P ( NABNCD  )  =0 • 0 
20  CONTINUE 
25  CONTINUE 

DO  35  I  I  =1 » IR1 

DO  30  I  I  I  =1  *  IR1 

I  I  I  I  I  =  I  I  I  +  I R 1* ( I  I -1  ) 

C ( I  I  1 1  I  )  =0 •  0 
30  CONTINUE 
35  CONTINUE 

DO  ino  J=1»JS 
Jn=IR*( J-1+JS#{K-1 ) ) 

ORS JK  =  0 • 0 
DO  70  I =1 » I R 
I JK= I+JK 
UBS  I JK=OBS (  IJK) 

OBSLOG=LOGF (OBSIJK ) 

DO  60  NB= 1 ♦ JS 1 
SIGMA=-1.0 

IF(J.EQ.JS)  GO  T0  40 
IF(NB.NE.J)  GO  TO  60 
SIGMA* 1.0 
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40  DO  50  N  A=  1  *  I R  1 
RhO=-l • 0 

I F ( I • £ 0 • I R )  GO  TO  45 
IFINA.NE.I )  GO  TO  50 
RHO= 1 • 0 

45  N Ari  =  NA  + 1  R  i  *  i  No- 1  ) 

G ( NAd ) =G ( NAd ) +RhO*6 I GMA*G6SLOG 
50  CONTINUE 
60  CONTINUE 

OBSJK=OBSJK+OBSIJ< 

70  CONTINUE 

DO  00  I  1=1  * IRl 

I  I  JK  = I I+JK 

DO  80  I  I  I  =  1 » IRl 

I  I  UK*  I  I  I+JK 

I  I  I  I i  =  I  I  I+IR1*(  I  I -1  ) 

I  I  1  I  I J* I  I  I  I I  +  IR1 IR1*( J-l ) 

Dim  I  J  =  -0do  (  I  I  J<  )  *0bi>  (  HUM  /03SJK 
i  F  (  I  I  •  £Q  •  I  I  I  )  3lIIIiJ  =  nIIIII  J+05S  (  I  I  J  K. ) 

C ( I  I  I  I i  ) =C(  I  I  I i I ) +61  I  i  I  I J 
a ( 1 1  I  I  I J )*dl  I  I  I  I  J 
80  CONTINUE 
90  continue 

too  continue 

aRITE  output  tape  output  *  2002 
CALL  NMPRINT(G*IR1JS1*1) 

CALL  MATINV(C*IR1*C*0*D*IR1) 

DO  170  Nd*l*JSl 
DO  160  ND= 1  *  Ji  1 
DO  150  N A  =  1  *  I R  1 
N  A  a  =  N  A  +  I Rl* ( Nb-1  ) 

DO  140  N C  =  1 » IRl 
NCD  =  NC+IR1*(ND-1  ) 

NA6NCD  =  NAB+IR1JS1*(\ICD-1 ) 

DO  120  I  1  =  1 » IRl 

I  I NCND= I I  +  IR1*(NC-1  +  IS1*IND-1 )  ) 

DO  110  1  I  I =1  * IRl 

NAI I IN6*NA  +  IR1*(  I  I I-1+IR1* (NB-1 )  ) 

I  I  I  I  I  =  I  I  I+IR1*<  I  1-1  ) 

P < NAdNCD ) =P ( NAdNCD )-ri(NAl I INa  )  *C (  I  I  I  1  I  )*6( IINCND) 

no  continue 

120  CONT INuE 

IF(Nd«NE»ND)  GO  TO  130 
NANCND=NA+ I  R  1  *  (  NC  -1  +  I  A  1*  ( ND- 1 ) ) 

P ( NABNCD) =P ( NABNCD ) +B(NAMC NO ) 

130  i*  (  NAdNCD  )  =0  (  NAdNCD  )  +  P(  NA0NCD  > 

PG(NAd)=PG(NAd)+P( NAdNCD ) *G ( NCD ) 

140  continuf 
i a o  continue 

1  40  CONT INUF 


II  /  -*»  / 
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170  CONTINUE 

DO  190  NA6  =  1 » IRUS1 
DO  180  NCD  =  1  * IR1JS1 
NA6NCD=NA6+IR1U31*(NCD-1 ) 

CHISQ=CHI5Q+G(NAB )*P(NA3NCD) *G(NCD) 

180  CONTINUE 
190  CONTINUE 

200  continue 

CALL  MATINV(Q,IR1JS1»Q»0*D*IR1JS1> 

DO  220»NAd=1 ♦ IR1JS1 
DO  210  NCD=1* IR1JS1 
NABNCD=NA6+IR1J31*(NCJ-1  ) 
Cnl3Q=CHlSQ-PG(NAB)*Q(NA6NCD)*PG(NCD) 

210  CONTINUE 
220  CONTINUE 

aRITe  OUTPUT  TAPE  OUTPUT, 2003 

WRITE  Output  TAPE  OUTPUT, 2004*NDEGREES»CHISG*IR*JS*<T 
300  END 
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oUBROuT  I  NE  T rtCB  YT  wO 

COM  i*  i  ON  /ALL/u6u(4U'j)  ,  I  JEN  T  (  10)  »  IR»US»K.T  »NDEGREES»  INPUT  » OUTPU  T 
J  I  MENS  I  On  0  (  16  )  »V  l  16  )  »6  (  16  )  »u  (  16  )  *b(  16  )  *  A  (  16  )  »D5T  AR  (  16  )  *VSTAR  (  16  )  ♦ 

*  GST  AR ( 16 ) »uSTAR ( 16 ) .BSTAR ( 1 6  )  .ASTAR ( 16 ) 

TYPE  INTEGER  OUTPUT 

2000  FORMAT ( 1H1»10A8) 

2001  FORMAT ( 26HOGOODMAN5  2X2XT  PROCEDURE  *////lH  ) 

2002  FORi*iAT(3HO  A»8X,3n  D  »10X,3H  V  *  1  oX  »3HX*X  »  1  uX  ♦  3H  Q  *]OX,3H  A  »1JX, 

*  3H4*Z  » 1  Ox  »  JH  G  *  1  OX  »  3H  U  ,10X*3HY*Y  »/lH  ) 

2003  FORMAT  (  I  3  »  9F 1 3  •  5  ) 

2  004  FORMA  T (  I  3  *  3 ( 2F13.6* 13X )  ) 

200^  F ORi,i A T  (///// 1  Vn  DcGREtS  OF  FREtDuM  » 1 1 X » 3 nX^X » 1 5X » 3 HZ *Z » 1 3X » 3H yu Y » 

*  15X.2H  R»8X*2H  o*8X,2H  T./lH  ) 

2006  FORMAT (6X,I3,6X,3F18.3»5X»3(8X,I2)) 

2007  FORMAT ( 3X ,3 ( 26X ,F13.5 ) ) 

a R I T E  OuTPuT  TAPE  0urpuT»20G0»IDENT 

ARITE  OUTPUT  TAPE  OUTPUT, 2001 

IF ( IR.Nc.2  )  GO  TO  100 

IFUS.NE.2)  GO  TO  100 

aRITE  OUTPUT  TAPE  OUTPUT, 2002 

SUM  1  =  SUM2  =  SuM3 -SuM4  =  SUMS  =  SUM6  =  SUM7  =  SUM  8  =  SUM  9  =  C . 0 
DO  20  K= 1 »<T 
KPART=IR*JS*(<-1  ) 

OBSl 1K  =  0BS ( 1+KPART  ) 

OBS21<  =  OBS( 2  +  KPART  ) 

0dS12<=06S( 3+KPART  ) 

OBS22K=OBS(4  +  <PART  ) 

D<=0BSllK*0oS?2K/06Sl2K/05S21K 

UK=1.0/OBS11K+1.0/OBS22<+1.0/08S1?<+1 .0/03S21K 

vk=dk*dk*uk 

A  <  =  1 , 0  /  V  < 

GK  =  LOGF ( OK  ) 

P<=1,0/UK 

BK= 1 .0/DK 

AK=3K*BK*UK 

HK= 1 • 0/AK 

UuM  i  =  ouM  1 +DX  D<*  »v  X 

SUM2=SuM2+DX*AK 

oUM3=5UM3+a< 

5UM4=SoM4+GK*GK*PK 
SUM5  =  SUV|5  +  GX*PX 
SUM6  =  SU’Y6  +  P< 

G  U  M  7  =  SUM  7 +6  X  *  B  K  *HX 
SUM8  =  SUM8+BX*HX 
SUM9=SUM9+HK 
DSTAR ( K ) =SUM2 / SUM 3 
VSTAR ( K ) =1 .0/SUM3 


Uk-Kk 


GSTAR ( < ) =SUM5/SUM6 
USTAR  (  0  =  1 .0/SUM6 
BSTAR(0=SUM8/SUM9 
ASTAR ( K )  =  1 .0/SUM9 
D  (  K ) =DK 
V  (  K ) =VK 
G  (  K ) =GK 
U(<)  =UK 
3  C  K  )  =  B  K 
A  (  K ) = AK 
20  CONTINUE 

SUMXSQ=SUMYSG=SUMZSQ=0.0 
<T  1  =  KT— 1 
DO  30  K  =  1 » K T  1 

XSQK  = ( D ( K+l ) -DSTAR ( O ) * ( D ( K+l ) -DSTAR ( < ) )/(V(<+l )+VSTAR(<) ) 

YSQK= (G( K+l ) -GSTAR ( K )>*(G(<+1>-GSTAR(0)/(U(K  +  1 )+USTAR (< ) ) 
ZSQK=(6(K+1 )-BSTaR<K) )*(B(K+1)-6STAR(0)/(A(K  +  1 )+ASTAR(K) ) 

WRITE  OUTPUT  TAPE  OU TPuT , 2 003 » < * D ( K ) » V ( O , XS OK ♦ B < O *  A ( < ) ,ZSQK » 

*  GU  )  *u  (  K  )  *  Y i>UK 
SUMXSQ=SUMXSG+XSGK 
SUMYSQ=SUMYaG+YSGK 
i>uMZSG  =  SUMZSG  +  ZbGK 

30  CONTINUE 

WRITE  OUTPUT  TAPE  OU  TPuT  .  2  004  ♦  K  T  »  D  <  XT  )  » V  (  KT  )  *  B  (  XT  )  » A  (  KT  )  ,G  (  <  T  )  * 

*  U  (  KT  ) 

WRITE  OUTPUT  TAPE  OUTPUT .2007. SUMXSutSUMZ5Q»SUMYSQ 
XXSQ  =  SUM1-SUM2*SUY,2/SUM3 
YYSQ*SUM4-SUM5*SUV|5/SwM6 
ZZSQ=SUM7-SUM8*SUM8/ UUM9 
WRITE  OUTPUT  TAPE  OUTPUT  *  2  OOA 

WRITE  OUTPUT  TAPE  OU TPUT » 2 006 *NDEGREES » XXSQ » ZZSQ . Y YSQ » I R ♦ JS » KT 
100  END 
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SUBROUTINE  kkk 

CGmm0N/ALL/UBo(40u  )  *  IOENT  (  10)  *  I  R  ♦  JS  »<T  »NDEGREES*INPUT  *OUTPUT 
COMMON/ SUMS /03SR ( 80 ) *OBSS ( 80 ) *08 ST ( 2 5 ) * 03 SR S ( 1 6 ) ♦ CPSTR ( 5 1 ♦ 

*  OBSST ( S )  .OBSRST 
TYPE  INTEGER  OUTPUT 

2000  FORMAT ( 1H1 » 10A8 ) 

2001  FORMAT  ( 39H0KUI_L6AC<  *  KUPPERMAN  *  AMD  <U  PROCEDURE  ) 

20  02  FORMAT  (/////  19rl  DEGREc.3  OF  FRtc.DOM  »  i  1  X  *  3  H2  *  I  ♦  1  7X  ♦  2  H  R*8X*2H  S»8X» 

*  2H  T  ) 

2003  FORMAT (/6X»I5»6X»F18.D*7X.3<6X.i2)> 
write  output  tape  output >2000* ident 
write  output  tape  output»200i 

CALL  GETSUMS 
TWOI =0.0 
DO  140  <  =  1  ♦  K  T 
OBSRSK=OBSRS( K  ) 

DO  130  J= 1 » JS 
OBSTRJ=OBSTR( J) 

JK  =  J+ JS* ( <— 1 ) 

0BSRJK=^6SR ( J<  ) 

DO  120  I =  1  *  I  R 
0BSSTI=0BSST( I ) 

I J=I+IR*( J-l ) 
uBSTl J=OBST ( I J) 

K  I  =K  +  K.T *  (  1-1  ) 

0BSS<I=0BSS(KI ) 

I  JK=I+IR*( J-1  +  JS*(<-1 )  ) 

OBSI  J<  =  OBS  (  I JK" ) 

TWOI sTWOI+OBSI JK*LOGF(OBSI JK*OBSSTI*OBSTRJ*OBSRS</OBSRST /OBST I j/ 

*  OBSSKI /OBSRJK ) 

120  CONTINUE 

i 3 0  CONTINUE 
140  CONT INUE 

TWOI =2.0*TWOI 

write  output  tape  output»2002 

WRITE  OUTPUT  TAPE  OU TPu T . 2 00 3  »NDEGRE ES * T wO I  *  I R * J S ♦ K T 
400  END 
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SUBROUTINE  INTRVALS 

COMMON/BCD/IBVECTOR(20) » JC VECTOR ( 2u)  .KDVECTO R(2u ) 

COMMON /ALL /OBS ( 40u ) ,IUENT( 10) » I R  * JS.KT  *NDEGREES* INPUT  »OUTPUT 

COMMON /FLAGS /OPT IOM 10) »LOGFLAG 

DIMENSION  0(6 ) »CHIT (2 ) *CHI 1 ( 2 ) »PHIW(2 ) »P (2 ) 

type  integer  output 

TYPE  INTEGER  OPTION 
DATA(P=0.95»0.99  ) 

2000  FORMAT ( 1H1 ♦ 10A8  ) 

2001  FORMAT ( 61H0 1 NTERACT I  ON  ESTIMATES  AND  CONFIDENCE  INTERVALS  IN  LOG  U 
♦NITS  ) 

2002  FORMAT (63H0INTERACTI0N  ESTIMATES  AND  CONFIDENCE  INTERVALS  IN  BASIC 
*  UNITS  ) 

2003-  FORMAT  ( ////  44H  B  C  D  ) 

2004  FORMAT (/9X»I2»l4X»I2»i4X»I2) 

2005  FORMA  I  l  /  /  /  /  34h  P  Chi  (  T  =  I  5  *43H  ) 


*  Chi  (  1 )  PHI U -  I  5  tlH) » / 1 H 

2006  FOkivmT  l9X»F4.2»3(9X*F15.4)  ) 

2007  FORMAT ( /////106H0  I  B  J  C  K  D  P 

♦TANDARD 

♦ERROR  CHI(T) 

*  PHI(W)  / * 1H  ) 

2008  F0RMAT(6I3»F9.2*2E15.S*4X»3(4X.2E10.2) ) 

2009  FORMAT ( /////106H0  I  B  J  C  K  D  P 

* 

*HI(T)  CHI  (  1  ) 

2010  FORMAT(6I3*F9.2*E15.3*15X»4X*3(4X*2E10.2) ) 
WRITE  OUTPUT  TAPE  OUTPUT  * 2 000 ♦ I  DENT 

IF (LOGFLAG.EQ.O  )  GO  TO  5 

kvR  I  Te  OUTPUT  TAPE  OOTPUT*200l 

GO  TO  10 

5  WRITE  OUTPUT  TAPE  0UTPUT*2002 
10  aRITE  OUTPUT  TAPE  OUTPUT»2003 
NW  =  0 

NPCD=n 

DO  15  1=1*20 

ib=ibvector ( ; ) 

IF(IB.EQ.O)  CO  TO  15 
JC  =  JCVECT OR (  I  ) 

KD  =  <DVcCTOR (  I  ) 

N6CD=N3CD+ 1 

NW  =N W+ ( IB-1 )♦( JC-1 ) * ( <D- 1 ) 

WRITE  OUTPUT  TAPE  OUTPUT , 2004 » I  9  ♦  JC »KD 
15  CONTINUE 

rtR I T E  OUTPUT  TAPE  OUTPUT »2u05 »NDEGREES*NW 
W  =  FLOAT F ( Nw  ) 

T  =  FLCATF( NDEGREES  ) 

DC  20  L  =  1  »  2 
Pl.  =  P(  I. ) 

CHITl  =  CHI  ( T .PL  ) 


) 


ESTIMATE  S 

...LIMITS...  /.P0X.77H 
CHI ( 1) 


ESTIMATE 

... LIMITS...  /.7?X**5HC 
PHI  (W)  / ♦  1  H  ) 
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CHI 1L  =  CHI  ( 1 .O.PL  > 

PHIWL  =  PHI  (  a1*  PL  ) 

wRITE  OUTPUT  TAPE  OU TPUT,2 006* PL »CH I TL*CH I  1L» PHI WL 
CHIT (L)=CHITl 
CHI  1 ( L ) =CH I iL 
PHI w ( L  )  =PHI wL 
20  CONT INUE 

IF(LOGFLAG.EO.O)  GO  TO  25 
WRITE  OUTPUT  TAPE  OUTPUT,2007 
GO  TO  30 

25  wRITE  OUTPUT  TAPE  OUTPUT»2009 
30  DO  100  ,N  =  1  *NBCD 
<D=<OVECTOP ( M ) 

<D  1  =  K.D-  1 
JC= JCVECTOR ( N ) 

JC1= JC-1 
I b= I BVECTOR ( N ) 

I  b 1  =  I  6- 1 

DO  70  <  =  1  * <D 1 

DC  60  J=l*JCi 

DO  50  I  =  1  * IB1 

1J<=I  +  IR*(J-1  +  JS*(K.-1  )  ) 

OBSI J*=ORS( I JK ) 

1dJCK=I3+IR*(JC-1+JS*( <-l ) ) 

ORRBC<=ORC ! IP JCK ) 

1BJK=IB  +  IR#(  J-1+JS*U-1  )  > 

0B56 JK=06S( I  3 J< ) 

IJCK=I+IR*(JC-1+JS*( <-l  )  ) 
uouICi<,=Obu(  IJCK  ) 

IJku  =  IflR#(J-l-t-Jb*(iCD-l)  ) 

ObSI JD  =  05S( I J<  D ) 

IBJC6D=lD+lR*(JC-i+Jo1<-(<D-l)  ) 

OBSriCD=OuS( I  3 JCKD  ) 

IdJKD= I8+IR*(  J-1+JS*(<D-1 ) ) 

0BS3 JD  =  035 (  IBJKO  ) 

I JCKD=I+IR*( JC-l+JS* ( KD-1 ) ) 

OBSI CD=OBS ( I JC  <0 ) 

aIBJCnD  =  LOGF( OBSI J<*0oS3CiC*0BSoJD*0bb ICO/C6SBJ</OBS IC</CrS I  JO/ 
*  CRSBCH) 

sj  I  o  J  C  4  =  1  •  u  /  U  ou  I  Ji^  +  1  •  u  /  uouoCA+i  «  u/UdS  I  C  <+ 1  •o/OBSBJK 
o  IoJCD=l«U/Oooi  Jji-i«0/U3SoCD-t-l»u/0DSICu+l  •^•/ubSBJD 
SIbJC<u=S0RTF(UIbJCK+JI3JCD) 

DO  45  L  =  1 ♦ 2 
PL  =  P ( L  ) 

w  (  1  )  =  AlbJCfvU-o  1  dJC6D*Ch  I  T  (  L  ) 

0(2 ) =2.0*AIdJC<D-Q(  1  ) 

0  (  3  )  =A  IRJCC-S  I  B  JCKD*CHI  1  (  l  ) 

Q(4)=2.0*AIt5jCKD-D(?  ) 

0(5)  =  A  I  B  J  C  K.  D  -  3  I  3JC<D*PHI  W(  L  ) 
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0(6) =2 • 0*A  I  dJCKO-Q ( 5  ) 

IF (LOGFLAG.tQ.O)  GO  TO  35 

•v  k  I  T  E  GuTPuT  TAPE  UuTPuT,2Gu8» I »IB»J*uC.6»6D.PL.AIBJCKD.SI8JC40»0 
GO  TO  45 

35  AIBJCKD=EXPF(AI5JCKD ) 

00  40  M  =  1  »  6 
0(M)*EXPF(Q(M) ) 

40  CONTINUE 

w K I T £  OUTPUT  Tape  OUTPUT »20lC* I » IB»J*JC *K»<0*PL »AIBJC<D *0 
45  CONTINUE 
50  CONT INUE 
60  CONT INUF 
70  CONTINUE 
100  CONTINUE 
ENO 
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1000 

1001 

1002 


400 


oudROUTlNE  RcADIN 

COMi»iOfc/ALL/OBS  (400  )  *  I  DENT  (  10  )  *  I R  ♦  JS  *KT  ♦NDFGREFS*  I  NPUT  ♦  OUTPUT 
COMMON/ BCD/ id VECTOR! 20) »JCVECT0R(2C )  ♦<DVECTOR ( 20 ) 


COMMON/FLAGO/OPT I  ON (  10) ♦LOGF LAG 
TYPE  INTEGER  OPTION 
TYPE  INTEGER  OUTPUT 
DIMENSION  INMAT (10) 

FORMAT ( 1 0A8 ) 

FORMAT (2014 ) 


FORMAT (1011) 

READ  INPUT  TAPE  I NPUT » 1000 ♦ I  DENT 
IF (IDE NT (1) •EG*ShTHAToALL)  GO  TO  400 
kc.hU  I  NPO  T  T  AP c.  I  NPu  I  »  1 0^2  *0P  T  I  ON 
READ  INPUT  TAPE  INPUT ♦1G01*IR*JS»<T 
READ  INPuT  T APE  I NPO T ♦ 1 000 » I NM A T 


i\R5T  =  IR->*-Jo'i'*GT 


READ  I  NP  J  T  T  A  Pc. 
READ  INPUT  TAPE 
RcAu  I NPu  T  T  APc 
READ  INPuT  TAPE 
N D E G R E F S  =  (  I  R-l ) * 


I  NP  J  T  ♦  I NMA  T» (CdS(N)  *  N  = 1 » NR  S  T ) 
INPUT  *  1001 ,IBVFCT0R 
INPUT* 1001 .JCVECTOR 
I  NPUT  » 1001 *KDVFCT0R 
(  JS-1 ) * ( <  T - 1  ) 


RETURN 

STOP 

END 
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subroutine  getsums 

COMmU(M/ALL/OB3(40j) *  I  DENT ( 10) * IR .JS.KT .NDEGREES *  I NPUT .OUTPUT 
COM^iON/SuMo/OBSRt  80  )  .GBSS(80  )  *OBST  (  25  )  .OBSRS  (  16  )  .OBSTR  (  5  )  » 

*  OBSST (5 ) .OBSRST 
OBSR  ST  =  Q • 0 
DO  30  K  =  1 *<T 
O8SRSK=0 •  0 
DO  20  J= 1 » JS 
OBSR JK  =  0  •  0 
DO  10  I  =1  ♦  I  R 
I  JK=I+IR*( J-1+J5*(<-1 )  ) 

UBS i JK=OBS ( I JK ) 

0BSRJK=08SRJ<+0BSIJK 
06SRSK=06SRSK+0BSI JK 
OBSRST =OBSRST+OBS IJK 
10  CONTINUE 

JK  =  J+ JS# ( K- 1 ) 

OBSR ( JK ) =OBSR JK 
20  continue 

OSSRS( < ) =OBSrSk 
30  CONT  INUE 

DO  60  J=1»JS 

OBST  RJ  =  0.0 

DO  30  1  =  1. IR 

OBSTl J  =  0.0 

DO  40  K=1.KT 

IJK=I+IR*(J-1+JS*(K-1 ) ) 

OBSI JK  =  03S( IJK  ) 

OBSTl J=0BSTI J+OBSI JK 
OBST RJ=OBSTRJ+GBS IJK 
40  CONT INUE 

I J=I+IR*( J-l  ) 

OBST ( I J ) =0B3T  I  J 
50  CONT INUE 

OBSTR( J)=OBSTRJ 
60  CONT INUE 

DO  90  1=1. IR 

OBSST 1=0.0 

DO  80  K= 1 ♦ K T 

OBSSKI =0.0 

DO  70  J=1 »JS 

IJK=I+IR*(J-1+JS*(K-1) ) 

OBSI JK=0BS( IJK ) 

OBSSKI =06SSK I +03S IJK 
OBSSTI=03SSTI+08SIJK 
70  CONTINUE 

K I=K  +  KT*( I -1  ) 

OBSS ( K I ) =OBSSK  I 
80  CONTINUE 

OBSST { I )=0BSSTI 
90  CONTINUE 
END 
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SUBROUTINE  NfVPRlNT  (  X  »N*M) 

COMMON  /  ALL/ods(4U  o)  *  I  uEN  T  (  10  )  ♦  IR»j5»X.T  *NDEGREES*  INPUT  »OUTPUT 
0  I  MENU  I  ON  X  (256)  *XNOR.4(2  5o  ) 

TYPE  INTEGER  OUTPUT 

2000  FORMAT ( /9X » 32HALL  NUMoERo  HAVE  riEtN  UIVIDED  dY  E11»1/*1H  ) 

2001  FORMAT (8X*16F7«3) 

NM=N*M 

XMAX=X ( 1 ) 

DO  10  K=2 ♦ NM 

IF(XMAX*XMAX.GT.X(<)*X(<) )  GO  TO  10 
XMAX=X ( K ) 

10  continue 

XMAXaAbSF ( XMAX ) 

XMAXLOG=LOGF (XMAX )/2. 302 5850930 
M  AX  L  OG  =  X  M  A  X  L  OG 

IF  (XMAX.LT.  1 .0  )  MAXLOo  =  MAXLOC--l 
XMAX=lO.O**KAXLOG 
DO  20  <= 1  *  NM 
XNORM ( K ) =  X ( K  )  /XMAX 
20  CONTINUE 

WRITE  OUTPUT  TAPE  OU TPUT » 2 000  •  XMA X 
DO  30  I  =  1  *  V 
J2  =N* I 
J 1  = J2-N+ 1 

WRITE  OUTPUT  TAPE  OU TPUT * 2 00 1  * ( XNORM l J ) » J  =  J 1 ♦ J2  ) 

30  CONTINUE 
END 


FUNCTION  PHI(W.PL) 

AREA=1.0+(PL-1.0)/2.0/w 

TEST=0. 

EP  =  .0000001 
IF ( AREA-. 5 )  3  *  4  *  5 
4  XORD=0. 

G0T016 

3  AREA1=1.-AREA 
C  =  -1. 

G0T07 

ARPA1  =ARFA 
C  =  1. 

7  X  =  ( AREA  1  -  .5)  *  2.5 

8  DIV  =  1.  +  .3275911  *  X 
E  =  1.0/  DIV 

6  =  ( ( ( ( .940646070*E) -1.287822453 )*E+ 1.2 59695 130) *E-.2 52 128668 )*E 

X+. 225836846 

oD  =  (((((  4.70323035  *  E  )  -  5.151289812  )*E  + 

X3. 77908539  )  *  E  -  .504257336  )  *  E  +  .22583o846  )  *  E 

AA  =  EXPF  (  X  **2  )  *  .88622692 
FX  =  (AREA1  -  1.  )  *  AA  *  DI V  #  2 •  +  S 
FXD  =  -  SD  *  .3275911  -  S  *  2.  *  X 

XNEW  =  X  -  FX  /  FXP 
I F ( ABSF  (XNEW-X)  -EP  )  6«  6*  14 

14  TEST-TEST+1. 

IF (TEST -40, ) 15*6*6 

15  X  a  XNEW 
GO  TO  8 

6  XORD  =  C  *  XNEW  *  1.414213562 

16  PHI =XORD 
END 
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function  chi ( t  ,pi j 
DIMENSION  X ( 1 3  »  2  ) 

iJATA(A  =  3*b4i*b*V>i»7«oi3*3i*4oo*ll*u7ufi*DDv*2»3253bi31*l.i‘,'iP« 

*  - u •  5349b7V 3? » “O •  1  ^3 0  i  179  7 *  J »v77V3Vbt> 7 ^*”0* lJw617d?6w#  1?24 ^ 5 i A? * 

*  6«635*1y*2lu*li#j43»l^#277»15«u36*l»uOj»3»2o946u7'r)f2*94^1b4^j* 

*  -0 •  2Vu3  lb:>y  4  *— 0*  d26:> *  0  •  9 1 1341798 1 “0 •  3423 1 56 13*0 *2v ^  1111 4®  ) 

L  =  1 

IF(PL»GT.0»97c  >  L  =  2 
I F ( T .GT • 5 • 5  )  GO  TO  ]0 
I  T  =  T 

CHI =X ( I T»L ) 

GO  To  30 
10  TRT  =  SORT f ( T  ) 

C  H  I  =  X  <  1  3  ♦  L  ) 

DO  20  K  =  1  »  7 
ChI=CHl/TRT+X(  1 3  — A  *  L ) 

20  CONTINUE 
CH I =T  *CH I 
30  CH I =SQRTF (Chi) 

END 
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T  ABl  F  1  -  OATA  FROM  KASTFN3AIJM  AMD  l.  AMPHIEAP  (1050)  -  2  X  3  X  5 

KASTFNRAUM-I  AMPH1FAR  PROCFOIIPE 


f.FU  PRSERVFO  FREQUENCY  FXPFC  TFD  FREQUENCY 


1 

1 

1 

SR. 000000 

54. 995074 

7 

1 

1 

75.000000 

78.004076 

1 

7 

1 

1 1 .000000 

17.331343 

7 

7 

1 

10.000000 

1 7. 618657 

1 

■> 

1 

5.000000 

6. 67363? 

7 

3 

1 

7.000000 

5.376363 

1 

1 

7 

40.000000 

48. 379034 

7 

1 

7 

5R.OOOOOO 

58. 620966 

1 

7 

7 

14.000000 

13.991513 

7 

7 

7 

1 7.000000 

1 7.008487 

1 

3 

7 

10.000000 

10.629453 

7 

3 

7 

R • 00^000 

7. 370547 

1 

1 

3 

33.000000 

34. 177160 

7 

1 

3 

45.000000 

43. 872840 

1 

7 

3 

ir.oooooo 

17.469760 

7 

7 

3 

77.O00O00 

77. 530740 

1 

3 

3 

15.000000 

14.403580 

7 

3 

3 

10.000000 

10.596420 

1 

1 

4 

15.000000 

17. 13208? 

7 

1 

4 

39.000000 

36. 86791 8 

1 

7 

4 

1  3.000000 

1 l. 079608 

7 

4 

77.000000 

73.920397 

1 

3 

4 

15.000000 

14. 78R31 1 

7 

3 

4 

1  3.000000 

18.211 6  R9 

1 

1 

5 

4.000000 

4.366700 

7 

1 

5 

5.000000 

4.633300 

1 

7 

5 

1  7 . ooo ooo 

1 3.07R277 

7 

7 

5 

15.000000 

13.971723 

1 

3 

5 

1 7.000000 

1 5. 555073 

7 

3 

5 

3.000000 

9. 444077 

1 TFRATIHN 


DFQRFFS  nF  FRFFnPM 


CH  I  -5QIJARr 


76 


3.  1  5RP45 


R  S  T 

2  3  5 
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’  ARI  F  2 


n  A  r  A  FPHM  KARTFNRA  M  amp  l  AMpHIFAR  (1R5Q) 


2  X  R  X  S 


PI  Af.KFTT  PR  nr  FOUR  F 
K=  1 

CONTRAST^ 

All  M'l^RFPS  HAVF  RFFN  HIVIDFO  PY 

?.  PRR  -1  . 

nTSPrPSTPN  MATRIX 

AM  MIIMRFRS  HAVF  RFFN  OIVIOFO  PY 

0.1 74  -n. 11 R 
-0.111  1.54  A 

K=  ? 

(  INTPASTR 

All  NIJMRFPS  HAVF  RF  F  M  HIVIDFO  PY 
0.7RR  -p.nqi 
niSPFPSION  MATRIX 

All  MtJMPrpe  HAVF  Rc  F v  niVIDFn  PY 
O.  1  SP  -0.  '"'9  i, 

-n.ooi  I.oap 

K=  R 

r—JTRARTR 

All  VIJmrf^c;  HAVF  Rrc»’  MVI'^'i  RV 

-0. 1  "M  -1 . R?7 

niSPFPSinN  MATRIX 

All  MIJmrcrs  HAVF  RFFN  DIVIOFO  PY 

1.RR5  -n.4PP 
—  0  •  4  p  5  P.?n? 


1 .0-00] 


]  ,04-non 


1  ,P-°01 


1  .0  +  r>no 


1  .r+oon 


1  .o-nni 
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TABLE  2  ( Continued) 


CONTRASTS 

ALI  NII^RFRS  HA VF  RFFN  OIVIOFO  RY  1.04-000 
-0.479  -1.117 

wspFRsiriN  matrix 

AM  NIIWRFPS  HAVF  RFFN  OIVTOF^  RY  l.o-oni 

?• 1 47  -0.001 
-0.001  7.004 

K  =  5 

CONTRASTS 

ALI  NUMRFPS  HAVF  RFFN  OIVTDFO  RY  l.O  +  '^OO 
0.000  —  1 . 9  54 
OT  SPFRSION  MATRIX 

All  NlJMRFRS  HAVF  RFFN  OIVIOFO  RY  1. 0  +  ^00 

0.6^0  0.000 

0.000  1.305 

7-VFf.TOR 

All  NIIMRFPS  HAVF  RFFN  OIVIOFO  RY  1.0  +  000 

-1  .  74 R  -5. 740 
F'-MATR  I  X 

ALI  NIJMR  FR  S  HAVF  RFFN  OIVIOFO  P,Y  1.0-001 

0.096  -O.ORR 
— 0 . OR  R  1.951 

SIJM(l  *L  ) 
9.530446 


OFGRFFS  OF  ^RFFDHM 


R 

2 


CHl-SOHARF 

3.  177641 

S  T 

3  5 
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TABIF  3 


DATA  FROM  KASTFNRAIJM  ANO  LAMPHIFAR  (1940) 


?  X  3  X  4 


Goodman  pRnr.rniiPF 
K=  1 

CONTRAST  s 

All  NIIMRFRS  hAVF  BEFM  OIVIDFO  BY 

P. 794  -7. IP? 

DISPERSION  MATRIX 

A U  NUMBERS  HAVF  B  F  F N  DIVIDED  BY 

3.774  3.479 

3.479  4.444 

K  =  7 

CONTRAST  S 

A I  I  NIIMRFRS  HAVE  RF  FN  DIVIDED  BY 

-3.91 R  -4. 1 73 

DISPERSION  MATP I  X 

All  NIIMRFRS  HAVF  BFFN  DIVIDED  BY 

7.474  7.740 

7.740  3.443 

K  =  3 

CONTRASTS 

A |  I  NIIMRFRS  HAVF  BEEN  DTVIDFD  BY 

-7.146  -6.06? 

DISPERSION  MATRIX 

Al  I  NIIMMTPS  HAVF  BFFN  OIVIDFO  BY 

7.197  1.667 

1.447  7.477 


1 .0-001 


1  .0-001 


1 .0-001 


1  .0-001 


I  .0-001 
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TABLE  3  ( Continued) 

K=  4 

CONTRASTS 

All  *\ll  I  MR  CR  S  HAVF  REFN  DlVlDFI)  BY  1,0-001 
-7.71?  -R.4BR 

dispersion  matrix 

All  NIJMRFRS  HAVF  BFEN  DIVIDED  BY  1.0-001 

7.145  1.7?? 

1.777  7.444 


K=  5 

CONTRASTS 

All  NUMRFRS  HAVF  RFFN  OIVIOFO  RY  1.0-001 
-0.749  -9.749 
DTSPFRSION  MATRIX 

AIL  NIJMRFRS  HAVF  RFFN  OIVTOFD  RY  1.0-001 

4  .VAR  1  .  ROB 
l.RRR  R.BRR 


G-VFCTOR 

All  NIJMRRR S  HAVF  RFFN  OIVTDFO  RY  1.0  +  000 
-7.00R  -4.51 1 
O-MATR  TX 

All  NHMRFPS  HAVF  RFFN  DIVIDED  RY  1.0-0"? 

5.4?4  3.RR4 

B.RR4  4.309 

SIJM  (  H4H  ) 
9. 5RR446 

R  S  T 

2  3  5 


DFGRFFS  OF  FRFFODM 


CHI-SOUAMF 


3.127441 
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TABIF  4 


DATA  FRO”  NORT ON1  (1945) 


7  X  2  X  12 


Goodmans  7X7xr  proofoijrf 


K 

D 

V 

X*X 

3 

A 

1 

0.05543 

0 . 0  0  0  R  4 

0. 29711 

IR. 039?7 

RR.597°4 

2 

0.03414 

0.00041 

0.07239 

77.65714 

742. 3R946 

3 

0.0471 7 

3.00066 

0. 71 299 

21.27105 

1 33.P7761 

4 

0.0670P 

0.00136 

0.96744 

16.10714 

°1 .40393 

5 

0.07715 

0.0O070 

0.O0713 

36.B3373 

375.06006 

4 

0.03947 

0.00045 

0.03231 

26.33333 

194.97354 

7 

0.03393 

0.00040 

0. 95997 

29.4734R 

300.12639 

R 

0.07353 

0.00131 

0.9799? 

13.60000 

44. P6733 

9 

0.07045 

0.00101 

0.  1  41 99 

14.1 9444 

40.97770 

i  r» 

0. 03314 

0.P0035 

1.69991 

30. 1 6000 

795.59509 

1 1 

0.01931 

0.00077 

7. 1 R260 

64.67600 

1 940. 99453 

1  7 

0.01 455 

0.00013 

7.451 5R 

6". 41 667 

1695.44771 

7*7 

G 

1 1 

V 

Y*Y 

0.77949 

-7. 99265 

0  .  ?  7  9  ?  3 

°. 30999 

0 . 00 1 R6 

—3.31 999 

°  .  3  1 6°  9 

n . OQ?  77 

0.16383 

-3.05499 

0. 7°  7  7  7 

0. 1 9R75 

0.76060 

-2.77976 

n . 3« 27  1 

°. 99999 

0. 1"675 

-3.40640 

0. ?76< 

0.02324 

0.71198 

-3.73212 

n. ?PP72 

0.  129,04 

1.00329 

-3.3835° 

0.34549 

1.16016 

0 .47970 

-7.61007 

0.2425? 

0.82690 

0.51606 

-2.45795 

0.20339 

0.41032 

0.67876 

-3.40652 

0.31397 

1.29140 

1.03793 

-4.00049 

0 . 6r'  046 

2.02544 

5.18475 

-4.10127 

0.46448 

7.37399 

OFGRFFS  OF  FRFFDOM 

X*X 

1*1 

Y*Y 

1  1 

7.45159 

5.19475 

7.373 

O  S  T 

2  2  12 
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TARI  F  5 


DATA  FROM  KASTFNBAIIM  AND  L  AMPHI  FAR  (1959) 


2X3X5 


T  NTFR  ATT  I  ON  FSTTMATFS  AND  OONFl  DFNCF  TNTFRVALS  IN  LOG  UNITS 


ROD 
7  3  5 


P 


rHKT=  fi ) 


r  h  i  ( i ) 


0.95 

0.99 


3.93R1 

4.4319 


1 .959R 
2.575R 


I 


R 

J 

r. 

K 

0 

P 

FSTIMATF 

ST  ANOAR  0 

ERROR 

7 

1 

3 

1 

5 

0.95 

1.056+000 

1.004+000 

7 

1 

3 

1 

5 

0.99 

1.056+000 

1.004+000 

7 

7 

3 

1 

5 

0.05 

7.663-001 

9.057-001 

7 

7 

3 

1 

5 

0.99 

7.663-001 

Q.057-001 

7 

I 

3 

7 

5 

0.95 

5.351-001 

9.46R-001 

7 

1 

1 

7 

5 

0.99 

S.R51-O01 

9.46R-00 1 

7 

7 

3 

7 

5 

0.95 

5.596-001 

R. 301-001 

7 

7 

3 

7 

5 

0.99 

5. 596-001 

R. 301-001 

7 

1 

3 

3 

5 

0.  95 

2.613-001 

9. 736-001 

7 

I 

3 

3 

5 

0.  99 

2.61 3-001 

9. 236-0O1 

7 

7 

3 

3 

5 

0.95 

3.703-001 

7.756-001 

7 

7 

3 

3 

5 

9.99 

3 . 7  0  R-  0  0  l 

7.756-001 

7 

1 

3 

4 

5 

0.  os 

2.037-001 

9.21 1-001 

7 

1 

3 

4 

5 

0.99 

2.037-001 

9.21 1-001 

7 

7 

3 

4 

5 

0.  95 

6.331-001 

7.605-0 U 

7 

7 

3 

4 

5 

0.99 

6 . 31 1-001 

7.605-001 
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TABLE  5  (Continued) 


PHI ( W®  81 

2.7344 

3.2272 


CHI(T) 

-2. 90*000  9.01*000 
-1.44*000  9.59*000 
-2.90*000  4.83*0"9 
-3.29*0"0  4 .83*000 
-9.14*000  4.91*000 
-9.66*000  4.49*000 
-2.71*000  9.99*000 
-9.16*000  4.28*000 
-9.98*000  9.90*000 
-9.89*000  4.40*000 
-2.68*000  9.42*000 
-3.U*9"P  9.85*000 
-9.42*000  9.89*000 
-9. 9?«-000  4.99*000 
-2.96*000  9.69*000 
-2.78*000  4.04*000 


...LIMITS... 

CHI ( II 

-o.H-OOl  02*000 
-1.59*000  9.64*000 

-1.01*000  2.44*000 

-1.57*000  3. 1 0*000 

-1.27*000  2.44*000 

-1.95*000  9.02*000 

-1. "7*000  2.14*000 

-1.58*000  2.70*00" 

-1.55*000  2.07*000 

-2.12*000  2.64*000 

-1.15*000  1.89*000 

-1.63*000  2  •  7*000 

-1.60*000  2. "1*00" 
-2.17*000  2.88*000 

-9.57-001  2.12*000 

-1.93*000  2.59*000 


PHKWI 

-1.69*000  9.40*00 
-2.14*000  4. **O*00 
-1.71*000  3.24*00 
-2.16*000  9.69*00 
-2.00*000  9.17*00 
-2.47*000  3.64*00 
-1.71*000  2.09*00 
-2.12*000  3. ?4*00 
-2.26*000  2.74*00 
-2.72*000  3.24*00 
-1.75*000  2.49*00 
-2.13*000  2.87*00 
-2.91*000  2.72*00 
-2.77*000  3.18*00 
-1.45*000  2.71*00 
-1.92*000  3.09*00 
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